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ABSTRACT 

The  compatibility  between  maneuver  load  control,  relaxed  static  sta¬ 
bility,  and  flying  qualities  requirements  is  investigated  in  this  report.  Three 
steps  were  involved  in  the  Investigation: 

1.  An  analysis  was  made  of  control  surface  combinations  and 

their  effectiveness  for  maneuver  load  control  when  used  with 
an  airplane  having  shortened  tail  length  and  reduced  tail 
surface  area. 


Control  system  configurations  were  syntliesized  that  minimize 
a  weighted  measure  of  change  in  drag,  wing  root  bending 
moment,  control  surface  activity  and  response  error  between 
a  Level  1  flying  qualities  model  and  the  actual  T-33  airplane. 


3.  A  direct  optimization  of  the  tail  length,  tail  area  and  control 
surface  deflections  required  to  obtain  a  compatible  com¬ 
promise  of  the  CCV  objectives  was  performed. 


The  results  show  that  reductions  in  maneuver  drag  and  wing  root 
bending  moment  can  be  achieved  if  sufficient  controllability  is  available  to 
generate  the  required  forces  and  moments  and  at  the  same  time  to  artificially 
compensate  for  the  lack  of  inherent  stability  of  the  vehicle. 
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The  research  documented  in  this  report  investigates  the  feasibility 
of  implementation  of  several  fundamental  Controlled  Configured  Vehicle  (CCV) 
concepts.  The  tail  length  and  area  of  a  typical  Air  Force  inventory  airplane 
(the  T>33)  were  decreased  to  reduce  the  static  stability  and  additional# 
active,  wing-mounted  control  surfaces  were  used  to  provide  a  measure  of 
control  of  the  wing  lift  distribution.  Control  system  design  concepts  were 
then  used  to  evaluate  the  effectiveness  of  the  maneuver  load  control  concepts 
and  to  produce  flight  control  system  designs  that  can  take  advantage  of  the 
additional  controllability  of  the  vehicle  and  at  the  same  time  provide  for 
adequate  flying  qualities.  There  is  very  little  doubt  that  the  geometry  and 
the  controllability  of  present  day  vehicles  can  be  significantly  improved  as 
soon  as  fly-by-wire  and  control  augmentation  are  fully  accepted  as  integral 
parts  of  a  flight  control  system. 

Once  feedback  and  command  augmentation  are  fully  accepted,  then 
the  full  potential  of  the  use  of  feedback  can  be  and  should  be  investigated  and 
developed,  for  this  knowledge  can  have  a  significant  impact  on  the  fundamen¬ 
tal  design  of  the  airframe.  For  instance,  inherent  longitudinal  or  directional 
stability  need  not  be  built  into  the  airframe  if  sufficient  control  surface 
effectiveness  and  power  exist  so  that  stability  can  be  maintained  by  feedback. 
This  is  generally  an  easy  requirement.  A  more  difficult  requirement  is 
that,  in  addition  to  stability,  the  vehicle  must  satisfy  flying  qualities  require¬ 
ments.  The  research  documented  in  this  report  places  special  emphasis 
on  the  flying  qualities  of  the  augmented  airplane. 

The  importance  of  the  use  of  active  controllers  in  addition  to  the 
conventional  moment  producing  devices,  i.  e. ,  elevator,  rudder  and  aileron, 
has  been  amply  demonstrated.  Operational  and  experimental  aircraft  now 
are  using  or  investigating  the  proper  use  of  active  X-force  control  devices 
(auto-throttle  systems),  Z-force  control  (direct  lift  flaps  or  spoilers)  and 
even  Y-force  control  (side  force  surfaces,  differential  throttle  or  differen¬ 
tial  drag  devices).  With  this  acceptance  of  additional  force  and  moment 
generating  devices,  it  becomes  important  to  investigate  their  use  in  the 
attainment  of  many  desirable  objectives  of  airplane  design,  such  as  control 
of  the  wing  lift  distribution  (to  reduce  the  critical  wing  root  bending  moment), 
structural  mode  control,  ride  qualities  improvement,  gust  alleviation,  and 
certainly,  the  dominant  factor  in  all  flight  control  system  design,  the  flying 
qvialities.  The  research  documented  in  this  report  emphasizes  the  use  of 
active  surfaces  such  as  a  direct  lift  flap  and  collectively  acting  ailerons  in 
addition  to  the  elevator,  to  produce  good  flying  qualities  in  addition  to 
reduced  wing  root  bending  moment  and  minimum  drag,  both  statically  (in 
trim)  and  dynamically,  during  maneuvering  flight. 

In  this  report,  the  T-33  airplane  was  used  as  the  object  of  the  study. 

A  fairly  modest  CCV  treatment  was  applied  to  the  vehicle:  the  tail  length 
and  size  were  assumed  variable;  the  flap  and  ailerons,  in  addition  to  the 
elevator,  were  assumed  to  be  actively  controlled  as  force  and  moment 


generating  devices.  The  object  was  generally  to  determine  the  tail  length, 
size  and  surface  motion:,  that  would  satisfy  the  generally  conflicting  require¬ 
ments  of  minimum  change  in  wing  root  bending  moment,  minimum  maneuver 
drag  and  minimum  deviation  from  Level  1  flying  qualities.  Much  more  could 
have  been,  and  eventually  should  be,  done.  The  usefulness  of  additional  force 
and  moment  devices,  such  as  canard  surfaces,  spoilers,  actively  controlled 
thrust,  sectioned  flap  and  sectioned  ailerons  could  have  been  included  in  the 
investigation.  They  were  not  included  for  three  reasons; 

1.  It  was  felt  that  the  configuration  selected  would  demonstrate 
the  feasibility.  Feasibility  rather  than  final  design  was  con¬ 
sidered  to  be  the  major  objective  of  the  program. 

2.  In  this  area,  which  represents  a  fairly  radical  departure 
from  conventional  design  and  development,  careful, 
measured  advances  should  be  made. 

3.  Time  and  money  for  an  extensive  investigation  were  not 
available,  and  in  fact  should  not  have  been  made  available 
at  so  preliminary  a  stage. 

A  second  area  in  which  the  study  was  restricted  was  the  area  of  more 
extensive  geometrical  changes  in  the  airframe  or  the  characteristics  of  the 
airframe.  The  tail  area  and  length  of  the  T-33  were  considered  changeable; 
these  changes  affected  the  stability  derivatives  of  the  airplane.  The  effects 
of  other  geometrical  parameters  were  not  investigated  but  eventually  they 
should  be  considered  in  the  CCV  context. 

Therefore,  the  study  was  restricted  to  feasibility  rather  than  a 
thorough  investigation  of  the  potential  of  CCV  concepts.  Yet  the  results  are 
both  promising  and  gratifying.  The  wing  root  bending  moment  can  be 
reduced;  the  drag  can  be  controlled,  and  the  T-33  airplane  can  be  made  to 
have  Level  1  flying  qualities  with  a  shorter  and  smaller  tail.  Therefore, 
from  the  CCV  point  of  view,  the  T-33  is  overdesigned;  weight,  drag  and 
structural  loads  can  be  reduced  by  using  a  feedback  augmented  flight  control 
system.  Since  the  demands  of  performance  right  now  result  in  dynamically 
less  well  behaved  airframes  requiring  augmentation  to  bring  their  dynamic 
behavior  up  to  acceptability,  this  report  takes  the  attitude  that  stability 
should  not  impose  strong  constraints  on  the  geometry  of  the  bare  airframe. 
Feedback  augmented  flight  control  systems  should  provide  the  required  sta¬ 
bility.  The  research  documented  in  this  report  shows  how  advanced  flight 
control  system  design  techniques  can  be  used  to  realize  many  of  the  CCV 
concepts. 

1.  2  APPROACH  AND  ORGANIZATION 

The  systematic  study  of  relaxed  static  stability,  maneuver  load  con¬ 
trol  and  attainment  of  good  flying  qualities  was  performed  in  four  major 
steps: 
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The  tail  size  and  length  were  varied  and  combinations  of 
elevator  motion,  inboard  flap  and  collectively  acting  ailerons 
were  used  to  force  the  vehicle  to  fly  "exactly"  as  a  flying 
qualities  model.  For  different  combinations  of  surface  usage, 
the  changes  in  drag,  wing  root  bending  moment  and  measures 
of  control  surface  effort  were  obtained  as  the  vehicle  maneu¬ 
vered  through  a  4  g  pullup  "exactly"  as  the  flying  qualities 
model.  These  calculations  verified  the  level  of  utility  of 
using  tlie  additional  surfaces.  This  research  is  described  in 
Section  II. 

A  feedback  flight  control  system  was  then  designed  to  minimize 
a  weighted  measure  of  incremental  wing  root  bending  moment, 
incremental  drag,  error  in  dynamical  response  between  the 
actual  aircraft  and  the  flying  qualities  model,  and  control 
activity.  Familiar  linear  optimal  control  techniques  were 
used  to  obtain  the  solutions.  The  effect  of  variations  of  the 
weighting  of  the  different  elements  in  the  performance  index 
was  investigated  to  some  extent.  The  purpose  was  to  verify 
that  the  different  weighting  could  be  successfully  "juggled" 
rather  than  obtaining  a  final,  or  "about-to-be-mechanized" 
design.  Section  III  of  the  report  describes  this  phase  of  the 
study. 

A  direct  parameter  minimization  of  the  performance  index 
was  tried.  The  purpose  was  to  try  to  directly  obtain  values  of 
tail  length,  size  and  control  deflections  that  would  minimize 
a  weighted  measure  of  incremental  wing  root  bending  moment, 
incremental  drag  and  minimum  deviation  from  Level  1  flying 
qualities  behavior  during  a  4  g  pullup.  This  effort  is  detailed 
in  Section  IV  of  the  report. 

Based  upon  the  more  than  twenty  years  of  experience  with 
feedback  fly-by-wire  systems  at  the  Cornell  Aeronautical 
Laboratory,  maximum  attainable  feedback  gains  for  the  present 
state  of  the  art  of  flight  control  system  design  were  estimated. 
Checks  were  made  to  assure  that  the  feedback  designs  defined 
in  Section  III  of  the  report  appeared  feasible. 
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SECTION  II 

ANALYSIS  OF  RELAXED  STATIC  STABILITY  AND 
MANEUVER  LOAD  CONTROL  SURFACE  REQUIREMENTS 


2. 1  INTRODUCTION 


The  first  logical  step  in  the  alteration  of  an  existing  airplane  to 
achieve  the  objectives  of  a  Control  Configured  Vehicle  is  the  determination 
of  the  effects  of  varying  the  tail  length  and  area  on  the  stability  of  the  vehicle 
and  the  ability  of  the  chosen  force  and  moment  generating  devices  to  maneu¬ 
ver  the  airplane  properly.  If  the  aerodynamic  characteristics  of  the  air¬ 
frame  do  not  inherently  provide  the  proper  forces  and  moments,  then  the 
surfaces  must  be  dynamically  moved  in  a  way  to  compensate  for  the  loss  in 
aerodynamically  generated  forces  and  moments.  The  analysis  of  geometry 
changes  and  surface  adequacy  was  accomplished  in  the  following  manner: 


1. 


The  equations  of  motion  of  the  T-33  airplane  as  a  function 
of  tail  length  and  area  were  developed. 


2. 


The  equations  of  motion  of  two  Level  1  aircraft  mathematical 
models  were  derived  from  MIL-F-8785(B). 


3. 


The  changes  in  drag  and  wing  root  bending  moment  were 
developed  as  a  function  of  tail  length,  tail  area,  and 
elevator,  flap  and  collectively  acting  aileron  deflections. 


4. 


Combinations  of  these  three  surfaces  were  then  deflected 
in  a  way  that  forced  the  T-33  to  respond  to  a  4  g  pullup 
command  as  the  flying  qualities  model  would  respond. 


5. 


Evaluations  were  then  made  of  the  effectiveness  and 
usefulness  of  the  surfaces  in  reducing  maneuver  drag, 
wing  loot  bending  moment  and  relative  surface  deflection. 


2.  2  DERIVATION  OF  T-33  EQUATIONS  OF  MOTION 
Equations  of  Motion  and  Dimensional  Data 


The  longitudinal  equations  of  motion  for  the  T-33  airplane  used  as  the 
example  in  this  CCV  study  are  as  follows: 


+ ^  Si  *  (^i.^  *  Kp  ^ej  j 


(1) 


P- 


i 
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Sj  -It  has  been  assumed  during  the  course  of  this  study  that  the  inboard 


flap  can  readily  be  changed  into  an  active  control  device,  requiring 
a  change  from  the  present  split  flap  arrangement  to  a  simple  flap 
with  the  following  characteristics: 


area,  each  side 
chord. 


Si-  =  15.  32  ft* 


s 


wing  station: 


=  1.8  ft 
0  ft  -  8.  5  ft 


.  It  has  also  been  assumed  that  the  T-33  ailerons  can  be  made  to  act 


collectively  to  a  longitudinal  stick  command.  These  outboard 
direct  lift  flaps  nave  the  following  characteristics: 


area,  each  side 

chord,  £g. 

wing  station: 


Sc.  =8.75  ft* 


=  1. 17  ft 

8.  5  ft  -  16.  0  ft 


Horizontal  tail:  area,  =  43.  5  ft* 
span,  b/f  =  15.  58  ft 
chord,  jSg  =  3. 12  ft 
taper  ratio  =  .  32 
42  =  5.6 


horizontal  tail  length  =  15.9  ft 


horizontal  tail  volume,  zsllL,:^  =  .  438 

H  s  z 


Sg  -  elevator 
area,  each  side 
chord. 


=8.7  ft* 


% 


=  .71  ft 


vertical  tail; 


Sy  =  22.  6  ft* 


area, 

vertical  tail  length =  16. 1  ft 

=  ,041 


Sv  — 


vertical  tail  volume  Vy  s 

^  h 
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Stability  Derivatives  of  the  Basic  CCV  -  T-33 

All  of  the  nondimens ional  stability  derivatives!  except  the  S/  and 
control  derivatives,  for  the  basic  T-33  were  obtained  from  Reference  1. 
The  control  derivatives  for  the  inboard  and  outboard  flaps,  5^  and  , 
were  estimated  from  Air  Force  Datcom  methods  (Reference  2)  using  the 
dimensional  data  above.  These  stability  derivatives  are  presented  in 
Table  I. 
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TABLE  I 

T-33  STABILITY  AND  CONTROL  DERIVATIVES 


Parameter 

FC-1 

FC-2 

Parameter 

FC-1 

FC-2 

0.019 

0.019 

0.  15 

0.  15 

0.  117 

0.  106 

6.50 

5.90 

0.0111 

0.0111 

0.362 

0.  343 

0.0176 

0.0173 

1,64 

1.64 

0.0408 

0.0408 

1.065 

1.065 

0.0703 

0.0703 

-0.010 

-0.009 

—  ^  4 

0.0192 

0.0192 

-0.690 

-0.590 

0.0384 

0,0384 

-3.30 

-3.  10 

0.06 

0.06 

-7,50 

-6.90 

0.019 

0.031 

-0.94 

-0.90 

C'Lt 

0. 142 

0.59 

-0.524 

-0.524 

0.00 

0.00 

< 

-0. 199 

-0. 199 

0.00 

0.00 

45.00 

10.00 

^1>AV 

0.00 

0.00 

Stability  Derivatives  in  Terms  of  and 


The  horizontal  tail  area  and  tail  length  of  the  T-33  were  altered  to 
Investigate  the  various  CCV  configurations.  The  following  parameters  are 
functions  of  the  tail  size: 


<2% 

K  1 


l: 


i>  t 


i 


i 
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Expressions  for  these  parameters  were  developed  in  terms  of  the  horizontal 
tail  area  ratio  and  tail  length  ratio  of  the  particular  CCV  configuration  under 
study  ?nd  the  basic  T-33  values,  where; 


=  — 
^  5, 


'CCV 


8A9IC  r-33 


■^Uccy 

r- 


tail  area  ratio 


tail  area  ratio 


(5) 


(b) 
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Methods  from  Reference  2  were  used  to  estimate  the  tail  contributions 
to  the  above  derivatives.  It  was  assumed  that  the  dynamic  pressure  ratio 
at  the  tail  was  1,0  as  the  horizontal  tail  is  more  than  five  feet  above  the  wing 
and  should  be  out  of  the  wake  of  the  wing.  However,  some  of  the  derivatives 
were  a  function  of  the  downwash  at  the  tail  which  was  estimated  from 
Datcom  methods  as: 


ll-  =  .484(.;j-^ 


(Datcom  Sec.  4.  4. 1) 


The  expressions  for  the  parameters  influenced  by  the  tail  size  are  listed 
below  with  the  corresponding  Datcom  sections  which  were  used  in  their 
derivation; 


Parameter  Datcom  Section 


'Ifel 


-w* 


Cm 


4.  15 

just  used  area 
and  length  ratio 
multiples 
of  basic 
derivatives 

4.5. 1. 1 

4.5.2. 1 

7. 4. 4.  2 

7.4.  1.2 


FC-1 

.0117 1:0015 
.out 
.3(,ZK/^ 

kt, 

3.ZS- 10.65  Kj^{Ky‘ 


FC-2 


.OiPi-.OOtS 
.01111/^ 
.3^3  K/j 


S.4S*[.819- 

US-'4.7Sk'^(kiy'^ 

Z.96'9.ssk:^((^j,y 
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Characteristics  oi  the  T-33  as  a  Function  of  and  A^ 

The  short  period  characteristics  of  the  T-33  were  calculated  for 
various  and  values,  to  see  what  the  bare  airframe  characteristics 
were  without  augmentation.  The  tail  ratios,  and  As^  ,  were  varied 
together  and  kept  at  identical  values.  A  simultaneous  reduction  in  both  tail 
length  and  tail  size,  with  corresponding  reduction  in  elevator  area,  could  be 
considered  to  be  a  fairly  severe  loss  of  stability  as  well  as  control  capability 
of  the  airplane.  However,  it  was  felt  that  tlie  loss  in  control  effectiveness 
and  control  power  could  be  maintained  with  an  all-movable  horizontal  surface. 
In  addition,  it  was  found  that  the  net  effect  of  inboard  flap  deflections  and 
outboard  flap  deflections  would  generally  be  in  a  direction  that  supplements 
the  pitching  moment  capability  of  the  elevator  during  pullup  maneuvers.  The 
best  values  of  A^4  and  investigated  separately  are,  of  course,  an  important 
CeV  concept,  but  for  this  study,  the  simultaneous  variation  of  these  param¬ 
eters  vividly  demonstrated  the  problems  and  design  principles.  At  these 
values  of  interest  they  both  had  approximately  the  same  effect  on  the  moment 
derivatives  which  directly  affect  the  short  period  mode. 

The  short  period  frequency  and  damping  and  for  these  configura¬ 
tions  were  estimated  by  the  following  approximations  (Reference  3); 


The  change  in  mass  and  moment  of  inertia,  ,  for  the  CeV  configuration 
with  /<4  and  Asi.  of  .  5  from  the  normal  T-33  were  estimated  at  a  negative 
three  percent  and  negative  ten  percent  respectively.  It  is  assumed  that  these 
vary  approximately  linearly  for  intermediate  values. 

The  short  period  characteristics  and  for  the  various  configura¬ 
tions  are  listed  in  Table  II  and  shown  in  Figure  1. 
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Derivation  of  W  ing  Root  Bending  Moment  and  Drag  During  Maneuvers 


To  determine  the  wing  root  bending  moment  (WRBM)  developed  in  a 
maneuver  it  is  first  necessary  to  calculate  the  spanwiss  lift  distribution  over 
the  wing  and  then  find  the  equivalent  force  and  moment  arm  to  which  this 
distribution  is  equivalent.  (It  was  assumed  that  the  lift  distribution  for  the 
wing  due  to  angle  of  attack  and  the  incremental  lift  distributions  due  to  flap 
deflections  are  linear  with  oc  ,  ,  and  So  and  are  simply  additive  with  no 

interference  effects  between  control  surfaces.  This  restricts  the  validity 
of  the  solutions  to  the  linear  range  of  {  let>j  £  15®)  and  |?/  6  25*.  This  is 
a  simplification  tliat  was  felt  to  be  justifiable  for  the  depth  of  investigation 
considered  in  this  study.) 


The  shapes  of  the  lift  distribution  curves  due  to  oc  ,  ^  ,  and  were 
estimated  from  methods  in  Reference  2  (Datcom  Section  6. 1.5.  1).  These 
are  shown  in  Figure  2.  The  area  under  each  curve  can  be  treated  as  the  total 
lift  increment  due  to  a  particular  ot  ,  ,  or  ^  deflection,  respectively.  The 

spanwise  lift  distribution  times  span  position  is  shown  in  Figure  3.  The  area 
under  each  curve  is  equivalent  to  the  incremental  WRBM  due  to  a  particular 
tL  t  Si  1  or  ^  deflection.  When  the  areas  under  the  curves  in  the  second  set 
are  divided  by  the  areas  in  the  first  set,  the  result  is  the  moment  arm,  or  the 
position  at  which  an  equivalent  force  will  produce  the  same  WRBM  as  the  lift 
distribution  did. 


Figure  1  ROOT  LOCUS  OF  CCV-T-33  SHORT  PERIOD  PCLESvs  TAIL  AREA  &  LENGTH  RATIOS 


Fi{;ii2.^2  WING  LIFT  DISTRIBUTIONS  FOR 
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The  above  procedure  was  carried  out  on  the  accompanying  figures 
with  a  planimeter  with  the  following  results: 


Lift  distribution  Lift  area  Lift  x  span 


due  to : 


position  area 
(VvRBM) 


Equivalent 
WRBM  arm 


.  173 


.244 


0788 


.414b  =  6.62  ft 
.323b  =  6.  16  ft 


.159 


.0874 


.548b  =  8.77  ft 


The  incremental  WRBM  derivatives  for  each  wing  are  given  below.  The 
moment  arm  associated  with  Mg  and  Mg^  are  the  same  because  the  lift  at 
06  =  0  has  the  same  distribution  as  the  lift  at  a:  51^  0. 


Mg  (zero  angle  of  attack)  ={kCi.^  ^  s)  *  (6.62  ft) 


=  *  (6.62  ft) 

=  (5.16  ft) 

=  *  (8.77  ft) 


FC-1 


FC-2 


41,970  ft-lb 
1,460,  500  ft-lb/rad 
357,  650  ft-lb/rad 
394,  750  ft-lb/rad 


12.  590  ft-lb 
397.  860  ft-lb/rad 
107,  300  ft-lb/rad 
118,  420  ft-lb/rad 


The  increment  in  WRBM  developed  in  a  maneuver  is  then  written  as 


OWRBM 


The  incremental  drag  developed  during  a  maneuver  is  defined  by  the 
following  nonlinear  equation: 


ADrag-^fsicn  fS^f^Cj,  j^oj 

'  ®  i  "u  -i  ^  fJ 


k^p  ICt,  a  +  a.  Si  ^  SoV. 
I 


I  iffducet^  drag  an  0ing) 

S  ) 

^  I  ^  Uivm  I'Tduced  drag  on  haritcnMf  HtU)  j 


In  some  instances  it  was  necessary  to  use  a  linear  dra^  equation  to  simplify 
the  calculations; 
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2.  3  DERIVATION  OF  AN  IDEAL  FLYING  QUALITIES  MODEL 

Good  flying  qualities,  as  defined  in  MIL-F-8785B,  is  one  of  the  major 
objectives  of  this  study.  Flight  control  system  designs  will  be  derived  that 
will  augment  the  airplane  in  such  a  way  that  the  augmented  vehicle  response 
will  be  identical  to  or  a  close  approximation  to  a  model  having  "ideal” 
flying  qualities.  In  this  section,  the  equations  of  motion  of  the  ideal  model, 
as  applied  to  this  particular  airframe,  the  geometrically  altered  T-33, 
are  derived. 

From  MIL-F-&785B  (Reference  4),  an  excellent  Level  1  longitudinal 
flying  qualities  airplane  would  possess  the  following  short  period  charac¬ 
teristics; 

.2 


Two  models  were  developed;  one  for  an  airplane  with  a  normal  lift  slope, 

Ci^  «  2fr  rad'l  (low  lift  model);  and  one  for  an  exceptionally  high  lift  air¬ 
plane,  Ci^  «  rad"^  (high  lift  model).  The  drag  polar,  velocity  derivatives, 
and  control  derivatives  were  assumed  to  be  the  same  as  the  normal  T-33. 
There  are  then  only  four  derivatives  that  need  to  be  found; 


It  was  first  assumed  that  =  Z  for  this  is  true  for  the  basic 
T-33  and  is  a  very  good  approximation  for  many  other  aircraft.  Then  the 
following  equations  can  be  written; 

^sp  ~  (ID 


~ 


-{3  Mi 


Z  = 


<2.  . 

F rom  the  relations  ci)^p  ~  n^lth  and  =.  7  for  an  ideal  Level  1  short- 
period  model  and  the  above  equations,  all  of  the  derivatives  of  the  ideal 
flying  qualities  model  can  be  evaluated.  Speed  stability  was  not  investigated 
in  this  study.  However,  the  phugoid  roots  of  the  models,  which  are  approxi¬ 
mately  those  of  the  T-33,  were  evaluated  and  were  also  Level  1. 

TABLE  III 

STABILITY  DERIVATIVES  AND  OTHER  CHARACTERISTICS 
OF  TWO  LEVEL  1  MODELS 


F  light 
Condition 

,  1/rad 
,  g/rad 
,  rad/sec 

a  -1 

I  S6C 


,  sec"  Wd”^ 
,  sec"  ^ 

, rad/sec 


Low  Lift  Model 


High  Lift  Model 


FC-1 

FC-2 

FC-1 

FC-2 

6.28 

6.28 

14.32 

14.32 

44.  12 

10.60 

100.7 

24.2 

6.  64 

3.26 

10.03 

4.92 

.7 

.7 

.7 

.7 

-2.214 

-.8237 

-5.045 

-L878 

-33.661 

-8.549 

-70.42 

-17.93 

-4.723 

-2.490 

-6.00 

-3.34 

-2. 362 

-1.245 

-3.00 

-1.67 

.06 

.  10 

.06 

.09 

.  11 

.04 

.  12 

.04 

The  above  models  are  labeled  on  Figure  4  along  with  the  unaug¬ 
mented  CCV  -  T-33  configurations.  It  shows  how  these  configurations  com¬ 
pare  to  MIL-F-8785B  short  period  frequency  requirements  for  Category  A 
Flight  Phases  (air-to-air  combat,  weapon  delivery,  etc.).  It  can  be  seen 
that  the  normal  T-33  (  =  1)  is  close  to  an  ideal  model,  while  the 

flying  qualities  deteriorate  quickly  as  and  are  reduced  below  .  85. 

Figure  5  shows  digitally  computed  time  histories  of  both  the  low  lift 
and  high  lift  model  responses.  The  primary  effect  of  n,/ci^  can  be  seen  in  the 
transient  of  the  pitch  rate  response  and  in  the  steady  state  value  of  Ao«  . 

High  v^foc,  means  higher  which  shows  up  as  a  smaller  overshoot  in  the 
pitching  rate  response  to  a  command  elevator  input. 
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2.  4  DERIVATION  OF  SURFACE  DEFLECTIONS  REQUIRED  TO  MATCH 

THE  MODEL  RESPONSE 

The  flying  qualities  model  defines  the  equations  of  motion  of  a  Level  1 
airplane.  Evaluation  of  surfaces  or  other  force  and  moment  generating 
devices  is  accomplished  by  obtaining  the  deflections  required  to  force  the 
CCV  -  T-33  airplane  with  altered  tail  length  and  size  to  respond  "exactly"  as 
the  flying  qualities  model  responds.  Three  control  surfaces:  elevator, 
inboard  flap  and  collectively  acting  ailerons  were  used  in  various  combina¬ 
tions  to  force  the  altered  airplane  to  respond  as  the  model. 

The  flying  qualities  model  represents  three  degrees  of  freedom  of 
motion;  if  fewer  than  three  independent  force  and  moment  generating  devices 
are  used  for  control,  the  response  of  the  flying  qualities  model  cannot  be 
exactly  reproduced;  only  one  state  variable  (and  derivatives)  per  controller 
can,  in  general,  be  made  to  respond  as  the  flying  qualities  model  would 
respond.  The  calculations  for  the  control  motions  are  done  in  the  following 
way; 


The  small  perturbation  equations  of  motion  for  an  aircraft  can  be 
written  in  the  state-vector  form 


i  - 


(15) 


where 

is  a  vector  representation  of  the  state  variables  of  the 
vehicle.  For  this  aircraft,  the  state  vector  0,/lvJ 

has  been  chosen. 


is  a  vector  representing  the  control  variables  of  the 
vehicle;  in  this  case  the  control  vector  elements  are  defined 
from  elevator  deflection,  Sg  ,  flap  deflection,  >  and 
collectively  acting  ailerons,  -  [Sg ,  S- ,  S^j 

F  is  a  square  (4  x  4)  matrix  of  dimensional  stability  derivatives 
of  the  airplane.  A  coefficient  of  this  matrix,  when  mul¬ 
tiplied  by  a  state  variable  and  an  appropriate  inertia  or  mass, 
represents  a  moment  or  force  applied  to  the  vehicle  due  to 
the  configuration  of  the  vehicle. 

(»  is  a  matrix  of  dimensional  control  derivatives  of  the  vehicle. 
A  coefficient  of  this  matrix,  when  multiplied  by  a  control 
surface  deflection  and  an  appropriate  mass  or  inertia, 
represents  a  force  or  moment  applied  to  the  aircraft  by  a 
surface  deflection. 
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Figure  6  MODEL  FOLLOWING  SYSTEM 

} 

t 

‘  Although  the  computation  to  determine  the  control  system  deflections^ 

U-(i)  actually  Involves  the  definition  of  a  control  law.  and  Incidentally  a 
control  system  deslgni  the  design  Is  generally  not  a  good  one  and  would  not 
normally  be  mechanized.  The  feedback  gains  can  be  destabilizing  and 
because  the  subset  x>2  of  the  state  vector  Is  not  required  to  be  fed  back, 
the  system  design  would  be  sensitive  to  poorly  known  or  varying  stability 
derivatives.  However,  If  the  same  number  of  controllers  as  degrees  of 
freedom  of  motion  were  available,  then  the  model  following  design  of  Figure 
6  Is  practical  and  any  feedback,  presumably  the  feedback  that  leads  to  a  stable. 
Insensitive  system,  would  form  the  basis  for  a  model  follov/lng  flight  control 
j  system  design. 


2.  5  RESULTS  OF  THE  STUDY 

The  curves  of  Figures  7,  8,  and  9  and  Table  IV  summarize  the  re¬ 
sults  of  using  the  model-following  technique  to  obtain  the  control  surface 
deflections  required  to  perform  a  4  g  pullup  which  would  be  performed  by  an 
airplane  whose  equations  of  motion  were  exactly  the  same  as  the  low  lift 
flying  qualities  model  at  FC-1  of  Table  II.  Figure  7  shows  the  change  In  drag 
obtained  by  using  only  the  elevator,  the  elevator  and  the  flap,  and  finally  the 
elevator  plus  collectively  acting  ailerons  as  a  function  of  tall  length  and  slza 
ratios,  and  .  The  curve  shows  that  the  maneuver  drag  is  reduced  both 
by  shortening  the  tall  and  reducing  the  tall  size  and  by  using  both  elevator  and 
flap  or  ailerons  to  perform  the  maneuver.  Drag  can  be  reduced  by  a' out 
13%  using  two  surfaces,  and  It  can  be  reduced  by  approximately  8.  5%  by  cut¬ 
ting  the  tall  length  and  size  In  half.  A  maximum  reduction  of  30%  can  be 
obtained  by  cutting  the  tall  length  and  size  In  half  and  using  both  elevator  and 
ailerons  for  maneuvering  control.  To  do  this,  however,  would  not  be  prac¬ 
tical,  for  as  shewn  by  Figure  9,  the  maximum  control  surface  deflection  of 
the  elevator  would  be  tripled.  Since  the  maximum  elevator  deflection  of  the 
T-33  Is  limited  to  -25.0*,  the  4  g  pullup  would  represent  the  maximum  at 
FC-1,  thereby  limiting  the  maneuvering  capability  of  the  airplane. 

21 


Figure  7  EFFECT  OF  TAIL  LENGTH  AND  SIZE  ON  DRAG  DURING  4  g  PULLUP 


SURFACE  DEFLECTION  (RAD) 


RATIO  OF  TAIL  LENGTH,  K^.  AND  SIZE,  K^;  K^  =  K^ 


Figure  9  EFFECT  OF  TAIL  LENGTH  AND  SIZE  ON  MAXIMUM  SURFACE 
DEFLECTIONS  REQUIRED  FOR  4g  PULLUP 
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The  reduction  in  the  maximum  wing  root  bending  moment  by  using 
auxiliary  surfaces  is  equally  pronounced  as  compared  to  changes  in  drag. 
Compared  with  the  nominal  size  tail  length  and  area,  the  wing  root  bending 
moment  decreases  by  approximately  9.2%  by  the  use  of  flap  and  elevator  and  by 
approximately  15%  by  using  elevator  and  collectively  acting  aileron.  The 
change  in  wing  root  bending  moment  also  decreases  as  the  tail  length  and 
area  are  reduced  but  the  change,  as  shown  in  Figure  8,  is  not  as  dramatic 
as  the  drag  effect. 

These  improvements  in  drag  and  moment  changes  cannot  be  obtained 
without  cost.  The  cost  is  obviously  in  the  deflections  of  the  control  surface 
required  to  perform  the  maneuver.  The  maximum  deflections  of  the  eleva¬ 
tor,  the  flap  and  the  ailerons  are  shown  in  Figure  9.  As  the  tail  area  and 
length  are  cut  in  half,  the  elevator  deflections  required  to  perform  the  4  g 
pullup  are  tripled,  with  some  decrease  in  maximum  elevator  required  when 
either  the  flaps  or  the  ailerons  are  used.  Because  the  elevator  deflection  ’s 
reduced  rather  than  increased,  the  flaps  and  ailerons  produce  a  beneficial 
effect;  they  aid  rather  than  fight  the  elevator  in  the  generation  of  the  response. 

influence  of  the  Flying  Qualities  Model 

The  flying  qualities  model  that  dictated  the  shape  of  the  response 
during  the  4  g  pullup  has  a  strong,  probably  the  most  important,  effect  on 
the  results.  Table  IV  below  summarizes  a  few  of  the  more  important 
results  using  the  high  Vf,/oi  model,  where  =  100.7,  as  compared  to  the 
lower  {v^/o^=  44.  12)  model  of  the  previous  analysis. 

TABLE  IV 

4  g  PULLUP  WITH  HIGH  MODEL 


Aircraft  Surface 

Configuration  Configuration 


A  Drag  A  WRBM  Max.  Surface  Deflec.(rad« 


Base  T-33 


only 
*  $i 


Sg^  *  Sg 


^g*Si  *Sg 


(lb) 

(ft-lb) 

907 

129, 100 

-.229 

2900 

139,  400 

-.  157 

2740 

158, 100 

-.  137 

2490 

184, 600 

-.  113 

814 

129, 100 

-.  381 

2970 

142,  300 

-.  361 

2800 

163, 300 

-.278 

2260 

179,  900 

-.  182 

+  .217 
+  .  128 


+  .228 
+  .  132 


+.  128 
+  .315 


+  .132 
+  .315 


By  comparing  the  surface  deflections  required  of  the  high  and  low 
requirements  (both  modelsi  high  and  low  Yi/ol>  ,  satisfy  the  flying  qualities 
requirements  in  MIL-F-8785B  equally  well)  there  is  a  significant  difference. 
For  the  low  vfoL  models  of  Figures  7,  8,  and  9i  the  flap  and  aileron  deflec¬ 
tions  were  negativei  which  decreat^ed  the  lift  on  the  wing  and  sigi'ificantly 
reduced  the  Wing  Root  Bending  Moment.  For  the  high  lift  model  at  FC-1 
whose  characteristics  are  tabulated  in  Table  III,  the  n/oc  of  the  model  was 
significantly  higher  than  that  of  the  T-33.  This  means  that  the  airplane 
performs  the  4  g  pullup  with  less  pitching  motion  than  it  normally  would  have. 
In  order  to  do  this,  the  effective  slope  of  the  lift  curve  must  be  increased. 

This  requirement  demanded  a  positive  (downward)  deflection  of  either  the 
flap  or  the  collectively  acting  ailerons.  An  increase  in  wing  lift  produces  an 
increase  in  wing  root  bending  moment  as  vividly  demonstrated  in  Table  IV. 

If  elevator  alone  is  used  to  generate  the  required  change  in  lift,  a  larger 
maximum  elevator  deflection  must  be  used  to  perform  the  4  g  pullup,  for 
l/'^2  larger  in  absolute  value  with  the  high  rt/oc  model.  When  two  control 
surfaces  are  used  for  maneuvering,  both  the  angle  of  attack  and  the  pitch 
rate  response  of  the  model  are  exactly  reproduced  by  the  T-33.  This  demands 
a  different  pitching -heaving  behavior  than  is  normally  obtained  in  a  T-33, 
causes  the  positive  flap  and/or  aileron  deflections,  and  increases  the  wing 
root  bending  moment.  Thus  it  can  be  concluded  that  relaxed  static  stability 
and  maneuver  load  control  are  incompatible  for  high  lift  models, 

Aa  the  ailerons  and  flaps  are  used,  the  deflection  action  requii “Aments 
of  the  elevator  are  significantly  reduced  and  the  elevator  deflections  bec^-'me 
less  strongly  a  function  of  tail  length  and  size.  This  is  to  be  expected  be¬ 
cause  heaving  (direct  lift)  requirements  are  predominant;  pitching  motions 
are  not  as  important  in  producing  lift  changes  and  therefore  become  a  less 
strong  function  of  tail  length  and  size. 

The  flying  qioalities  then,  have  a  very  pronounced  effect  on  the  maneu¬ 
ver  load  alleviation  requirements  of  a  Control  Configured  Vehicle.  The 
requirements  to  be  satisfied  by  the  CCV  should  be  as  low  as  possible  con¬ 
sistent  with  available  elevator  power  if  wing  trailing  edge  lift  modulating  sur¬ 
faces  are  used.  This  will  generally  require  negative  (T.E.  up)  lifting  control 
surface  deflections  that  will  aid  the  elevator  and  at  the  same  time  redistribute 
the  lift  on  the  wing  in  a  desirable  way  to  reduce  the  wing  root  bending  moment. 
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SECTION  III 

THE  SYNTHESIS  PROBLEM 


3.  1  INTRODUCTION 

In  the  previous  sectioni  the  T-33  airplane  was  forced  to  respond  as 
the  flying  qualities  model  responds.  The  tail  length  and  size  were  varied  to 
determine  the  effect  of  these  parameters  on  the  drag>  wing  root  bending 
moment  and  on  the  surface  deflections.  Several  combinations  of  candidate 
surface  configurations  inherent  to  the  T-33  airplane  that  can  be  considered 
for  maneuvering  were  considered.  These  surfaces  were  the  elevator,  the 
flaps  (assumed  to  be  a  simple  flap  capable  of  both  positive  and  negative  de¬ 
flections),  and  collectively  acting  ailerons.  Other  surfaces,  such  as  the 
addition  of  spoilers  or  canard  surfaces  were  briefly  considered  then  dis¬ 
carded,  because  it  is  felt  that  these  additional  force  and  moment  generating 
devices  should  not  be  added  to  the  airframe  unless  the  existing  surfaces  prove 
incapable  of  doing  the  job. 


% 


In  the  previous  section,  the  approach  taken  was  one  of  analysis. 
Parameters  were  changed  and  the  effects  were  noted.  In  this  section  the 
problem  of  synthesis  is  investigated.  Since  flying  qualities  requirements  are 
fairly  broad  and  many  different  configurations  will  yield  a  Level  1  aircraftr 
the  requirement  that  the  airplane  respond  exactly  as  the  model  is  much  too 
stringent  a  requirement.  More  emphasis  should  be  put  on  drag  and  wing  root 
bending  moment  reduction  for  this  was  felt  to  be  actually  more  important  than 
satisfying  the  ultimate  in  flying  qualities. 


To  satisfy  the  conflicting  requirements  of  minimum  drag  and  wing  root 
bending  moment,  a  performance  index  was  formulated  that  included  quadratic 
measures  of  the  error  in  dynamic  behavior  between  the  actual  aircraft  and  the 
flying  qualities  model,  the  maneuver  drag,  the  change  in  wing  root  bending 
moment  and  the  control  surface  motions.  The  quadratic  performance  index 
is  an  indirect,  rather  than  direct  measure  of  the  design  objectives.  The 
modeling  error,  the  drag  and  the  wing  root  bending  moments  are  minimized 
relative  to  each  other  in  a  way  that  produces  a  most  useful  kind  of  solution  to 
the  problem.  The  control  motions  are  relatively  smooth  and  well  behaved  and 
the  control  effort  and  maximum  deflections  are  managed  by  the  judicious 
choice  of  weighting  parameters  in  the  performance  index.  The  resulting 
control  law  is  linear  for  a  linearized  description  of  the  airplane  dynamics, 
and  these  f  )ntrol  laws  are  of  the  type  most  likely  to  be  actually  mechanized 
on  an  aircraft.  The  parameters  of  the  closed  loop  system  will  then  yield 
results  that  will  indicate  those  stability  derivatives,  such  as  and 
that  a  CCV  aircraft  might  inherently  possess.  The  performance  index  is  of 
the  general  form 
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where 


L  -  matrix  of  dimensional  stability  derivatives  of  the  flying 
qualities  model 


AD-Bxi-hIu  =  linearized  expression  for  the  change  in  drag 


AV12BM  =  t^u  =  linearized  expression  for  the  change  in  wing  root 
bending  moment 


M 


^  =  control  deflections;  =  [^e  •  >  ^e>  ] 

i!'  =  state  vector  of  the  airplane;  'X,^-  ^AV,  A9,  ef. ,  AcoJ 

(/% 


The  solution  is  constrained  by  the  equations  of  motion  of  the  T-33 
airplane. 


X  =  -h  du. 


(24) 


where  P”  and  9  are  the  matrices  of  dimensional  stability  and  control  deriva¬ 
tives  of  the  T-33  airplane  given  in  Table  I. 


The  solution  to  the  problem  posed  above  will  yield  the  motions  of  the 
three  control  surfaces  as  a  function  of  the  state  vector  <4  -  -  that  will 
minimize  the  performance  index.  The  matrices  Q  ,  ^  ,  T  and  B  express 
relative  emphasis  placed  on  the  requirements  to  minimize  model  following 
errors,  drag  changes,  wing  root  bending  moment  changes  and  control  motions. 
These  weighting  parameters  are  adjusted  relative  to  each  other,  the  absolute 
numbers  are  comparatively  meaningless  except  in  very  simple  cases. 


The  problem  of  minimizing  a  quadratic  performance  index  is  a  well 
established  method  of  flight  control  system  synthesis  and  literally  hundreds 
of  papers  and  reports  have  been  written  on  the  subject  since  the  technique 
was  formulated  and  popularized  by  R.E.  Kalman  (Reference  7)  and  S.S.  Chang 
(Reference  8).  Later  reports,  like  Reference  9,  established  the  relationships 
that  exist  between  the  performance  index  form  of  solution  (called  linear 
optimal  control)  and  the  more  conventional  control  system  synthesis  tech¬ 
niques,  like  root  locus  methods.  Reference  9  gives  many  examples  of  real¬ 
istic  flight  control  applications  and  the  theory  is  very  briefly  summarized. 
The  solution  to  the  linear  optimal  control  problem  yields  the  following  con¬ 
trol  lav/  which  minimizes  the  performance  index  of  Equation  (23) 


(25) 


where  P  is  the  positive  definite  symmetric  solution  to  the  matrix  Riccati 
equation 


0  =  PF^F  -  PGR'^C,*P  +  Q  . 


(26) 
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As  shown  in  Appendix  I  this  equation  satisfies  the  general 
Hamiltonian  equations 


-G,2''g' 


^(o)  = 
UO)  «  to 


that  are  necessary  and  sufficient  for  the  solution  of  this  problem.  The 
mathematics  of  this  problem  will  not  be  discussed;  however,  a  brief  intro¬ 
duction  to  one  of  the  more  important  closed  form  solution  techniques  is  given 
in  Appendix  I.  A  few  characteristics  of  the  solution  that  will  result  when  this 
technique  is  applied  are  discussed  below. 


Guaranteed  Stability  -  If  the  bare  airframe  is  unstable,  the 
choice  of  positive  definite  weighting  matrices  Q  ,  V  ,  T 
and  Q  will  always  yield  stable  linear  solutions  to  linear 
problems.  This  does  not  mean  that  the  resultant  closed  loop 
system  will  necessarily  be  stable.  If  an  error  between  a 
model  and  the  actual  airplane  is  minimized,  the  closed  loop 
aircraft  response  can  be  unstable  if  the  model  is  unstable, 
but  the  error  between  the  aircraft  and  the  model  will 
approach  zero  asymptotically. 

The  solution  will  generally  yield  a  closed  loop  system  that 
has  a  smooth  and  well  behaved  response  in  the  state  variables 
(or  errors)  included  in  the  performance  index.  These  states 
will  generally  respond  more  quickly  than  the  open-loop 
aircraft  and  exhibit  little  overshoot  to  an  initial  condition  or 
command  input. 

The  solutions  as  a  function  of  the  weighting  matrices  exhibit 
no  surprises.  A  series  of  solutions,  which  require  at  most 
a  few  moments  of  digital  computation  time,  quickly  establishes 
the  trends  of  the  solutions  as  a  function  of  the  v/eighting 
matrices.  Engineering  judgment  based  upon  the  knowledge 
of  the  limits,  capabilities  and  flying  qualities  of  the  airframe 
is  used  to  adjust  the  weighting  matrices  to  rapidly  arrive  at 
acceptable  solutions.  The  actual  numerical  value  associated 
with  the  performance  index  is  a  very  poor  substitute  for 
knowledge  of  the  airframe  stability  and  control  and  flying 
qualities  requirements. 


Input  Design 

The  linear  optimal  control  problem  described  above  to  obtain  the 
minimum  integral  of  drag  change,  wing  root  bending  moment  change  and 
minimum  integral  error  squared  during  a  transient  or  set  of  initial  con¬ 
ditions  describes  only  the  feedback  or  regulator  part  of  the  solution  to  the 
problem.  The  input  or  command  gains  must  also  be  defined.  The  problem 
could  have  been  formulated  as  a  model  following  problem  as  sketched  below: 


with  a  performance  index 
oo 

0  ' 


(32) 


but  the  regulator  or  closed  loop  part  of  the  system  is  not  influenced  by  the 
model;  the  model  is  an  uncontrollable  part  of  the  system  and  appears  only  in 
the  feedforward  or  command  portion  of  the  system.  It  was  felt  to  be  more 
realistic  in  terms  of  an  operational  system  to  include  the  model  in  the  perfor¬ 
mance  index  as  a  restraint  on  the  feedback  or  regulator  part  of  the  system 
and  separately  compute  the  feedforward  gains  to  yield  a  good  quasi-steady 
state  match  of  the  model,  at  a  time  t  -  Z  seconds  after  the  applied  command. 
At  this  time  the  short  period  had  responded  but  there  was  no  significant  speed 
change.  This  was  felt  to  be  realistic  in  terms  of  the  majority  of  maneuvering 
requirements  of  existing  fighter  aircraft. 


Three  control  surfaces  are  used  to  control  the  three  degrees  of 
freedom  of  motion  of  the  vehicle  so  the  problem  can  be  exactly  solved.  The 
command  input  gains  are  obtained  by  solving  for  the  values  of  the  control 
vector  Uf.  in  the  following  equation 
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“c  =  KsT'e'l^i„(t^)-[F-ciK],i„(t,) 

(vhire-  ii  -  Z  sec  a-fiar  4he  applied  pilot  input- 


This  guarantees  that  the  states  of  the  augmented  airplane  have  the 
same  values  (  q  ,  u  >  V  ,  0  ,  q  ,a;,V)as  the  model  2  seconds  after  the 
command  input.  The  trajectory  of  the  response  between  ^  =  0  and  t  -  Z  sec 
will  be  different  and  this  difference  will  be  a  function  of  how  closely  the  re¬ 
quirements  of  minimum  dynamic  error  between  the  model  and  the  actual  air¬ 
craft  were  met  compared  to  the  other  minimization  requirements  of  the  per¬ 
formance  index.  An  alternate  way  to  consider  the  command  part  of  the 
system  is  to  connect  the  stick  only  to  the  elevator,  but  this,  as  will  be  shown 
later,  produces  a  less  desirable  solution. 

Example  Solution 

A  fairly  wide  range  of  linear  optimal  control  solutions  was  run  to 
investigate  the  feasibility  of  obtaining  less  drag  and  wing  root  bending  mo¬ 
ment  during  dynamic  maneuvers.  The  important  considerations  were  integral 
of  drag  or  absolute  magnitude,  and  peak  wing  root  bending  moment.  Additional 
evaluations  were  based  on  maximum  control  deflections  and  control  action, 
(integral  of  the  square  of  the  control  deflections)  as  well  as  evaluations  of  the 
time  histories  of  the  response  of  the  system. 

Solutions  were  easily  obtained  that  yielded  good  results  and  could  also 
be  mechanized  without  difficulty.  An  example  is  given  by  the  performance 
index 


W  ^  SS  5/  f  5/ 1  dt 


This  performance  index  weights  each  portion  of  the  flying  qualities 
error,  the  drag  changes  and  the  wing  root  bending  moment  changes  in  the 
same  order  of  magnitude  and  is  a  very  straightforward  way  to  select  the 
weights  in  the  performance  index.  The  solution  produces  a  feedback  gain 
matrix 


(dV) 

(&9) 

(e) 

9.9  X  10'^ 

9.6  X  fo'* 

~.2S 

-1.0Z 

.124- 

-.196 

4,4x  /O'* 

'.in 

5.S& 

which  indicates  that  feedback  from  pitch  rate  and  angle  of  attack  changes  to 
the  three  control  surfaces  are  the  only  significant  feedback  requirements. 

The  gains  are  quite  reasonable  and  relatively  easy  to  mechanize.  The  closed 
loop  system  matrix  becomes 


F-flA?  = 


•32.17 

+3.82  X  {O'* 

32.17' 

0 

0 

1.0 

0 

2.637*10'^ 

*2.01  X  10'^ 

-6.823 

-76.S9 

*6.19x10'^ 

.9976 

-109_ 

The  feedforward  gains  obtained  from  Equation  (34)  are  shown  in  the 
sketch  below  relative  to  the  basic  T-33,  which  requires  an  elevator  deflection 
of  Sg  =  -  .  158  rad  to  obtain  a  quasi-static  change  in  normal  acceleration  of 
4  g  at  the  flight  condition  under  study,  P  C-1. 


Figure  10  COMMAND  INPUT  MECHANIZATION 


Figure  11  shows  the  responses  of  the  flying  qualities  model  and  the  optimal 
system  to  a  4  g  command  input.  The  response  of  the  augmented  vehicle  is 
similar  to  that  of  the  model  and  would  most  likely  be  considered  to  have 
Level  1  flying  qualities.  The  motions  of  the  control  surfaces,  although  ini¬ 
tially  abrupt  because  no  actuator  dynamics  were  included  in  the  simulation, 
are  not  considered  excessive. 


Table  V  shows  the  results  of  three  different  solutions  with  the  basic 
T-33  airplane  while  Table  VI  shows  the  feedforward,  and  feedback  gains 
and  closed  loop  eigenvalues  for  the  three  systems.  The  tables  show  that  the 
maneuvering  drag  and  the  wing  root  bending  moment  can  be  reduced  by  a 
reasonable  amount  through  the  reduction  of  the  tail  length  and  area,  yet  the 
control  system  required  to  still  give  good  Level  1  flying  qualities  is  not 
overly  complex  and  can  be  mechanized  without  difficulty.  The  only  potential 
problem  is  that  about  4  times  as  much  elevator  deflection  is  required  to  per¬ 
form  the  4  g  pullup  when  the  tail  length  and  size  are  cut  in  half.  This  may 
be  solved  by  replacing  the  elevator  with  an  all-movable  horizontal  surface. 
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TABLE  V 

COMPARISON  OF  CCV  PERFORMANCE  FOR  4  g  PULLUP 
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io‘’ 

108 
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Weighting  1 
Parameters 

% 

0 

10*' 

3 

0 

0 

V 

0 

10 

0 

0 

T 

0 

1 

1 

10 

0 
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Drag  (lb) 

2200 

2400 
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1830 
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112, 500 

127, 600 

127,  200 

-.0758 

-.085 

.265 

.262 

— 

.059 

-.078 

-.075 

z 

— 

-.  147 

.019 

.017 

0 

--- 

.  127 

.403 

.311 
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Effect  of  Command  Input  Mechanization 

The  linear  optimal  control  optimization  study  of  this  section  produced 
feedback  gains  but  did  not  specify  the  command  or  feedforward  gains  and  this 
produces  somewhat  of  a  problem. 

One  obvious  way  to  specify  command  gains  is  to  design  the  flight  con> 
trol  system  as  a  model  following  system.  A  second  way  is  to  calculate  the 
feedforward  interconnecting  gains  such  that  the  resulting  control  effectiveness 
matrix  would  be  the  same  as  that  of  the  model.  A  third  way  would  be  to 
command  only  the  elevator,  and  altering  the  feedforward  gain  or  the  "gear 
ratio"  in  such  a  way  that  the  augmented  aircraft  maintains  a  quasi-static 
4  g  response  to  the  same  stick  deflection  as  the  basic  T-33  airplane.  A 
fourth  way  of  designing  the  command  input  portion  of  the  system  is  given  by 
Equation  (34),  where  the  gains  are  calculated  such  that  the  augmented  air¬ 
craft  response  has  the  same  state  values  as  the  model  2  seconds  after  the 
command  input  is  applied  by  the  pilot. 

In  this  section,  a  comparison  is  made  of  two  command  input  designs, 
command  to  the  elevator  alone  and  command  to  produce  a  state  match  at  f  - 
2  seconds.  The  feedback  used  was  generated  by  the  third  solution  to  the 
linear  optimal  control  problem  of  Table  V,  in  which  the  weighting  param¬ 
eters  were  =  10^,  t  =  10“2,  R  =10°I. 

Comparisons  of  the  responses  of  the  systems  are  shown  in  Figure 
12.  The  most  significant  difference  between  the  two  system  responses  is 
in  the  control  deflection  time  histories.  The  design  that  commanded  only 
the  elevator  input  requires  significantly  larger  peak  and  steady  state  control 
deflections  to  obtain  the  4  g  pullup  with  significantly  larger  increase  in  drag 
during  the  maneuver.  The  deflections  required  of  the  elevator  command 
system,  however,  are  in  the  right  direction,  inboard  flap  down,  outboard 
flap  up,  to  produce  significantly  less  wing  root  bending  moment  change 
during  the  maneuver. 

There  are  actually  a  near- infinite  number  of  ways  that  the  control 
surfaces  can  be  connected  to  the  stick  command  input.  At  least  four  loeical 
ways  are  mentioned  above,  and  each  has  its  advantages  and  disadvantages. 

Yet  all  .with^e  possible  exception  of  the  model  following  arrangement)  can 
considered  to  be  solutions  to  the  linear  optimal  control  problem  solved  in 
this  section,  for  any  initial  set  of  control  deflections  can  be  thought  of  as  a 
set  of  initial  conditions  of  the  state  vector,  since  the  control  law  u.  ^ 
directly  relates  the  control  deflections  and  the  state  vector. 
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RESPONSE  WITH  STICK  COMMANDS 
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Figure  12  COMPARISON  OF  RESPONSE  WITH  TWO  DIFFERENT  INPUT  DESIGNS 
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Since  CCV  concepts  involve  reduced  tail  lengths  and  size,  using  addi¬ 
tional  control  surfaces  to  aid  the  resulting  less  effective  elevator,  the  inputs, 
i.e.,  the  deflections  of  the  various  surfaces  required  to  perform  the  maneuver^ 
are  critical  to  the  CCV  design  concept.  These  deflections  along  with  the  tail 
length  and  area  ratios  are  directly  considered  and  optimized  to  achieve  the 
desired  objective  of  good  model  following  with  minimum  change  in  drag  and 
wing  root  bending  moment. 


3.  3  LATERAL-DIRECTIONAL  RESULTS 


Introduction 


Although  most  of  this  study  dej.ls  with  a  CCV  for  the  longitudinal 
degrees  of  freedom  of  motion,  it  is  desirable  to  check  similar  vertical  tail 
variations  on  the  lateral-directional  behavior  of  the  airplane.  A  brief  inves¬ 
tigation  was  conducted  similar  to  the  linear  optimal  control  solution  for 
longitudinal  CCV  described  earlier  in  this  section. 


The  lateral-directional  equations  of  motion  are: 

V*'  yr->' 

0  1  0 

0  L' 


-p 

• 

r 

N'p 


•r 


L 

0 

0  0 

4 

r 

1 

Sf 


(38) 


As  was  done  with  the  longitudinal  modes,  it  is  necessary  to  derive 
relationships  of  the  lateral-directional  derivatives  with  the  altered  vertical 
tail  areas  and  lengths  of  the  CCV  and  the  normal  T-33  derivatives. 


The  basic  T-33  derivatives  for  the  following  flight  condition  were 
obtained  from  Reference  1 

\4  =  805  ft/sec 

h  =  23,  000  ft 

9  =  372  Ib/ft* 
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=  2,  5  rad/sec 
=.25 
=  .  25  sec 
‘2's  =  oo 

0//5  =  .  9 

The  performance  index  included  quadratic  measures  of  the  error  on 
the  dynamic  behavior  between  the  CCV  and  the  model,  and  the  control  sur¬ 
face  motions:  oa 

0 

where  both  of  the  weighting  matrices  were  set  equal  to  identity  matrices.  A 
more  comprehensive  study  should  consider  tail  loads. 

Seven  CCV  configurations  with  various  values  of  K/^  and  were 
evaluated.  The  results  are  listed  in  Table  VII. 


Evaluation  of  Results 

All  of  the  modal  characteristics,  though  not  identical  to  the  "excellent" 
model,  are  still  Level  1  according  to  specifications  in  Reference  4.  However, 
there  are  other  factors  which  show  that  there  is  a  limit  to  the  amount  of 
vertical  tail  reduction  that  is  allowable. 

First  of  all,  ^  feedback  gains  much  greater  than  5.  are  not  realizable, 
so  configuration  7  with  a  realistic.  Also  the  amount  of 

rudder  deflection  to  hold  a  constant  sideslip  may  be  a  limiting  factor.  From 
MIL-F-8785B,  an  airplane  must  be  able  to  bold  an  approximate  15  degree 
sideslip  in  case  of  an  extreme  crosswind  landing.  The  amount  of  control 
deflection  necessary  to  trim  each  CCV  configuration  in  a  steady  level  15 
degree  sideslip  is  also  listed  in  the  following  table.  Again,  configuration  7 
is  unrealistic  as  it  would  call  for  a  57  degree  rudder  deflection. 

It  can  be  seen  from  configurations  2  and  4,  and  3  and  6  that  the  effects 
of  and  are  almost  identical  and  stability  derivatives  depend  on  the 
value  of  the  product  ’  /^  )•  The  limiting  value  for  the  CCV  -  T-33  for 
( ^  most  likely  0.  1.  F rom  the  table  it  can  be  seen  for  this  value 
(configurations  3  and  6)  the  gains  are  about  5.0  and  the  rudder  deflection 

to  hold  the  maximum  sideslip  is  approximately  20  degrees.  Though  the 
vinaugmented  configuration  6  is  statically  and  dynamically  unstable: 
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RESULTS  OF  LATERAL-DIRECTIONAL  STUDY 


Para. 

"Excel¬ 
lent 
Model 
(Ref.  6) 

Basic 

Unaug¬ 

mented 

T-33 

1 

Augmented  CCV 

2  3 

Configurations 

4  5 

6 

7 

-- 

-- 

1.0 

5. 

0.  1 

1. 

.5 

.2 

.  1 

-  . 

1.0 

1. 

0.  1 

.5 

.5 

.5 

.5 

2.5 

2.2 

2.45 

2.41 

1.90 

2.41 

2.28 

1.87 

1.50 

.25 

.09 

.  167 

.  171 

.546 

.  184 

.283 

.555 

.773 

'^e 

.25 

.3 

.250 

.250 

.248 

.250 

.249 

.249 

.249 

% 

oo 

150. 

685. 

882. 

138. 

905. 

430. 

130. 

66.4 

.90 

L6 

.92 

.92 

1.28 

.92 

.99 

1.29 

1.78 

Response 

F  eedback 
Gains 

sjfi 

.356 

.252 

.  167 

.534 

.325 

.202 

.159 

-- 

-- 

0 

0 

0 

0 

0 

0 

-.001 

V(= 

-- 

-- 

-.051 

-.053 

-.053 

-.051 

-.053 

-.053 

-.054 

^cjr 

-- 

-- 

-.008 

-.005 

0 

-.024 

-.021 

-.013 

.003 

SrhS 

-- 

-- 

-.044 

.964 

5.800 

.980 

2.671 

5.727 

8.956 

-- 

-- 

0 

.001 

.082 

.001 

.012 

.080 

.201 

^rh 

-- 

-- 

.010 

.010 

-.014 

.013 

.015 

-.020 

.031 

Sfir 

-- 

-  - 

-.03V 

-.  100 

-2.447 

-.119 

-.485 

-2. 385 

-5.548 

Sf^  for 
Steady  15* 
Sideslip 

-- 

-12.4 

-12.4 

-9.7 

-9.3 

-10.4 

-10.3 

-9.6 

-9.5 

$^.  for 
Steady  15* 

0.93 

0,93 

9.5 

-19.8 

9.3 

-1.0 

-20.2 

-57.3 
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(the  two  real  Dutch  roll  roots  are  -1.  7  and  +1. 05  and 
Tg  =  .34  sec,  =  -  11.4  sec) 

with  a  simple  feedback  system,  this  configuration  can  be  made  to  fly  with 
Level  1  handling  qualities.  Therefore  a  corresponding  and  possibly  neces¬ 
sary  reduction  in  vertical  tail  area  and  length  is  compatible  with  reductions 
in  the  horizontal  tail.  Also  with  a  smaller  tail,  the  vertical  tall  loads  will 
be  reduced,  resulting  in  a  possible  reduction  in  structural  stiffening  and  a 
further  weight  saving. 


SECTION  IV 


DIRECT  OPTIMIZATION 

4.  1  INTRODUCTION 

In  Section  II  of  this  repo.-ti  several  control  surface  configurations 
were  investigated  in  order  to  obtain  measures  of  change  in  dragi  wing  root 
bending  moment  and  surface  deflection  as  the  T-33  airplane,  with  altered 
tail  length  and  area,  was  forced  to  respond  "exactly"  as  the  flying  qualities 
model  would  respond.  This  study  revealed  the  effectiveness  of  surfaces 
other  than  the  elevator  in  aiding  the  elevator  to  perform  the  pullup  maneuver. 
Then,  in  Section  III,  the  control  system  v/as  obtained  that  would  minimize  a 
measure  of  the  change  in  drag,  wing  root  bending  moment  and  error  between 
the  actual  aircraft  and  the  model.  Two  values  of  tail  length  and  area  were 
used  to  show  the  effect  of  this  parameter  on  the  resulting  solutions  but  no 
direct  effort  was  made  to  optimize  tail  length  and  area.  It  was  also  shown 
that  the  input  design  has  a  great  effect  on  the  resulting  system  behavior. 

In  this  section,  the  geometry-dependent  characteristics  a^e  treated 
directly.  The  objective  is  to  determine  optimum  tail  length,  tail  area  and 
surface  deflections  that  would  minimize  the  trim  drag  and  wing  root  bending 
moment  and  the  change  in  drag,  wing  bending  moment  and  the  model- response 
error  for  a  4  g  pullup.  The  problem  is  open  loop  in  the  sense  that  no  feed¬ 
back  will  be  directly  obtained;  instead,  the  geometrical  parameters  will  be 
optimized,  and  the  control  deflections  required  to  minimize  drag  and  wing 
root  bending  moment  in  trim  and  in  maneuvering  flight  will  be  treated. 

It  was  felt  that  the  dynamic  optimization  would  provide  the  optimum 
tail  length  and  area  and  once  vhese  were  obtained  the  deflections  of  the  sur¬ 
faces  could  be  calculated  to  maintain  minimum  drag  and  wing  root  beading 
moment  in  trim.  Then  the  required  command  inputs  obtained  from  the 
dynamic  optimization  would  complete  the  design  requirements.  In  actual 
practice,  the  two  parte.,  static  and  dynamic  optimization  were  done  separately 
yet  concurrently,  so  the  static  optimization  did  not  use  the  values  of  and 
1:2  obtained  during  the  dynamic  optimization  part  of  the  study.  It  would  be 
a  relatively  simple  matter  co  repeat  the  static  optimization  design  for  any 
value  of  Kn  and  , 

4.  2  STATIC  TRIM  OPTIMIZATION 
Formtilation  of  the  Problem 

It  is  desired  to  minimize  drag  and  WRBM  while  trimming  the  CCV  with 
the  elevator,  inboard  flap  and  outboard  flap.  The  following  functional,  F  , 
is  formed  containing  a  measure  of  drag,  WRBM,  and  Lagrange  multipliers 
times  the  three  longitudinal  trim  equations  (  level,  1  g  flight  with  constant 
thrust  and  velocity); 


F  =•  *  (P*'<ig)  ■*■  kij^  {\/^R.B^)  A/  (Trim  EquatioHi) 

uJhem  Vn^  -  Op^^  S;  *  Op^  %  tCp^^  % 

■^  K  '  Ci,*  Si  'Cp  Spf  ^ 

«)  ^ 


(41) 
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WRBM  s:  Mg  4- Mg  ai  +Mg  Si  S 

6  A  °Si  *Q 

(43) 

Trim  Et^aaHans 

h  (';r = s  'V '  ^ 

(44) 

(45) 

^gfo  =  6^  4  C>mJ^  *  ^j»jr  4j  *" ^Wg 

\  0  c6  *o  h  ' 

(46) 

Taking  derivatives  of  F  with  respect  to  the  unknown  variables  and 
setting  them  equal  to  zero,  we  obtain  seven  equations  v/hich  can  be  solved 
simultaneously  for  a  unique  solution  for  each  set  of  weighting  constants, /(i' 
and  ,  on  the  drag  and  WRBM< 


(51) 


The  seven  unknowns  are:  the  trim  a  t  3i  >  §o  •  ^e  •  Lagrange 

multipliers  X,  ,  X2  •  values  of  the  latter  three  parameters 

have  no  significance  as  the  equations  of  motion  are  just  added  to  insure  that 
the  aircraft  is  trimmed.  Only  the  relative  difference  in  the  weighting 
constants,  and  ,  is  significant,  and  if  both  are  increased  by  the  same 
multiple,  the  resulting  solution  will  be  the  same  except  that  the  Z/  will  chang 
by  the  same  multiple.  The  actual  values  of  and  have  no  meaning  other 
than  they  do  cause  different  minimum  drag  and  WRBM  solutions  to  occur. 

Results  and  Evaluation 


The  only  CCV  -  T-33  configuration  investigated  was  the  one  with  the 
normal  T-33  tail  size,  =  1.  Other  coi^figurations  would  involve 

similar  results.  The  significant  solutions  obtained  are  summarized  in  the 
following  table: 

TABLE  VIII 

SUMMARY  OF  TRIM  OPTIMIZATION  RESULTS 


OL 

(rad) 

(rad) 

(rad) 

(rad) 

Trim  (Level 

1  g  Flight ) 
(Ib)^ 

Trim 
WHIM 
(ft- lb) 

lO.xlO"^ 

1. 

.  168 

-.426 

2,069 

-4.071 

16, 313 

-580,498 

5.xl0"^ 

1. 

.020 

...093 

.380 

-.683 

4,  278 

-52,  575 

4.4x10""^ 

1. 

.007 

-.063 

.230 

-.381 

3,218 

-16,  312 

4.xl0"'^ 

1. 

-.003 

-.041 

.  117 

-.  154 

2,499 

18,  316 

3.6x10’^ 

1. 

-.015 

-.015 

-.013 

.  107 

1,953 

58,  221 

l.xlO”^ 

1. 

-.079 

.  130 

-.760 

1.589 

6,  956 

284,  551 

.  IxiO'"^ 

1. 

-.  110 

.200 

-1.  103 

2.294 

9,374 

392,  622 

1. 

1. 

-.  113 

.207 

-1.142 

2.373 

9,  645 

404, 669 

trim  with  normal 
T-33  (no  S;  or 

Sj,  available) 

-.001 

-.010 

— 

r  •  • 

1,723 

40, 998 
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It  can  be  seen  that  the  trim  WRBM  can  be  eliminated  completely  with 
inboard  and  outboard  wing  flap  deflections.  However,  this  is  done  at  the 
expense  of  increased  trim  drag.  In  fact,  it  appears  as  though  all  of  the  trim 
solutions  with  and  develop  more  drag  than  the  basic  T-33  solution.  Also 
the  trim  flap  deflections  are  relatively  high  for  the  reduced  WRBM  solution. 
For  example,  the  solution  which  results  in  an  approximate  50%  reduction 
in  WRBM  (  kT,  =  4.xl0"'^),  requires  almost  10  degrees  of  flap  deflection.  This 
would  severely  limit  the  amount  of  flap  deflection  left  for  maneuvering.  Each 
CCV  under  study  must  be  evaluated  separately  to  determine,  for  its  particular 
mission,  whether  the  reduced  trim  WRBM  and  the  resulting  decreased 
structural  weight  is  worth  the  increased  trim  drag  and  control  deflections. 


4.  3  DYNAMIC  OPTIMIZATION 

Derivation  of  Necessary  Conditions 

The  longitudinal  small  perturbation  equations  of  motion  being  con¬ 
sidered  are; 


«  =  >  ^(o)  “  0 

y  -  A  jb  f  Bic 

'ii  *  {v,e  tcY  -  perturbed 
fi  ^  $  KlV  ■  horizontal 

a  =  (Sg  ,  -  control  ps 

9  ,  -  measure  ( 


-  perturbed  state  vector 

-  horizontal  tail  geometrical  parameter 

-  control  perturbation  vector 

-  measure  of  state  at  final  time 


It  is  desired  to  find  a  step  perturbation  in  the  control  vector  H  -  hi  and  the 
tail  parameters  ,  such  that  during  the  maneuver,  which  transforms  the 
CCV  from  the  initial  state  to  some  final  state  y  (if)  at  a  prescribed  final  time, 
if  ,  the  motions  of  the  airplane  will  be  close  to  an  ideal  flying  qualities 
model  and  the  drag  and  WRBM  will  be  minimized.  In  the  following  discussions 
we  are  concerned  only  with  the  drag  and  WRBM  developed  during  the  maneu¬ 
ver  so  we  will  use  just  the  linearized  incremental  expressions: 


A  Vraq  -  AP  =  ^  +  (f^)  u 

«  Aw  -  f  rnltL 


Also  required  in  the  above  expressions  are  the  sensitivities:  and  ^'^Isbj 

which  are  the  solution  of: 
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The  above  expressions  are  reduced  somewhat  in  the  calculations  because 
many  of  the  matrices  contain  many  zero  locations.  The  function  dependence 
of  the  stability  derivatives  were  reduced  to  linear  relations  in  and  . 
This  gave  the  following  matrices  for  the  CCV  -  T-33  at  Flight  Condition  1: 
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Letting  the  parameter  vector  be  defined  as: 

P  =  (^^  .  ,  Si  ,  ,  a., ,  ,  %y) 

the  conjugate  gradient  algorithm  can  be  written  as  follows: 

A 

(i)  Guess  an  initial  value  for  P  and  call  this  P, 

(ii)  Compute  J  and  Vp  7 C^i)  using  the  previous  expressions 
and  let 

(iii)  Determine  the  constant  multiplier  tt^  of  the  corrections 
by  a  one  dimensional  search: 

(iv)  Define  the  new  parameter  vector:  •‘oOiO’^ 

(v)  Evaluate  the  updated  J  and  test  to  see  if  it  has  converged 

to  a  minimum  value.  If  it  has,  a  solution  has  been  obtained. 
If  it  has  not  converged,  proceed  to  the  next  step 

(vi)  Compute  the  gradient  Vp  j(P^) 
and  new  corrections: 

/A,  IVpJ(P/)t 

‘‘‘  '  ‘ 

(vii)  Repeat  steps  (iii)  through  (vii)  until  a  minimum  J  has  been 
reached. 


4.4  RESULTS 

Results  of  a  few  representative  runs  of  the  conjugate  gradient  program 
are  presented  in  Table  IX.  Two  different  initial  value  sets  were  run  with 
'/arious  weighting  constants: 
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Set  1:  initial  1.0 »  -  .S 

Set  2;  initial  .01,  kfi  =  .01 

The  maneuver  performed  was  a  4  g  pullup  in  2  seconds  with  9  (i^)  =  2  rad/sec, 
Htf)  =  0  rad/sec. 

As  can  be  seen  in  the  table  of  results,  the  solutions  are  not  too  pro¬ 
mising.  The  conjugate  gradient  method  as  formulated  here  yields  many 
local  minima.  In  fact,  for  each  different  set  of  initial  guesses  for  ^  , 

step  I  etc. ,  a  different  solution  for  the  minimum  of  the  performance 
ndex,  T ,  was  obtained.  Also  the  size  of  the  control  step  inputs  relative 
to  each  other  had  a  much  more  pronounced  effect  on  the  solution  that  was 
obtained  (the  initial  guess  in  each  case  was  Sl  =*15  rad  step)  than  to 
changes  in  tail  area  or  tail  length.  In  many  solutions  the  tail  parameters 
increased  in  size.  For  Set  2  all  of  the  solutions  remained  with  the  smallest 
tail  and  even  though  the  conditions  at  2  seconds  were  matched,  the  airplane 
was  still  unstable,  and  was  diverging  rapidly  at  that  time. 

However,  the  conjugate  gradient  method  may  still  be  of  some  use. 

One  remedy  to  the  non-uniqueness  problem  may  b  ‘  to  use  a  fewer  number  of 
unknown  variables.  By  reducing  the  number  of  control  step  inputs  as  vari- 
bles,  there  would  most  likely  be  fewer  local  minimurr.s  and  more  importantly, 
the  solution  would  have  to  be  reached  by  the  changing  of  the  tail  parameters, 
and  ^  rather  than  through  the  use  of  control  deflections. 

More  research  in  this  area  is  necessary  and  perhaps  a  different 
minimization  technique  less  prone  to  local  minima  should  be  investigated, 
such  as  quasilinearization  or  Kalman  filtering.  A  better  choice  of  cost 
function  may  also  be  found. 


52 


TABLE  X 

ATTAINABLE  FEEDBACK  GAINS 


Variable 

Surface 

Elevator 

Fi  .p 

Other  Wing 
Surface 

Rudder 

Aileron 

/av 

— 

— 

— 

— 

AoL 

rad/rad 

±5 

±5 

»  • 

±4 

• 

d 

rad/(rad/sec) 

:fe3 

±3 

±Z 

•  •• 

±3 

9 

rad/ rad 

±5 

±4 

±2 

... 

•  •  - 

rad/rad 

.... 

... 

±3 

±8 

-P 

rad/(rad/sec) 

— 

— 

— 

±2 

±3 

raa^raa 

— 

— 

±4 

r* 

rad/(rad/sec) 

— 

— 

±3 

— 
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’T'*??^??t$??%#C8WlV'W«KSWJV»wn*Hr'-«w^ 


5.1 


SECTION  V 

CONTROL  SYSTEM  MECHANIZATION 
INTRODUCTION 
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In  previous  sections,  the  need  for  adequate  control  power  was  stressed. 
If  the  tail  size  were  reduced,  including  the  elevator  size,  it  is  fundamental 
that  larger  elevator  surface  deflections  would  be  required  to  maneuver  the 
airplane.  Flap  and  collectively  acting  ailerons  help  the  elevator  produce  the 
required  pitching  mcixients  but  not  to  a  completely  compensating  amount,  so 
there  is  a  limit  to  the  extent  that  the  elevator  size  can  be  reduced.  Control 
surface  power  is  a  fundamental  limitation  associated  with  a  CCV. 

A  second  fundamental  limitation  associated  with  a  CCV  is  the  physical 
limitation  associated  with  feedback  control  mechanization.  Sensor  and 
amplifier  noise  and  structural  flexibility  limit  the  amount  of  feedback  that 
can  be  applied  to  an  airplane.  Table  X  summarizes  some  of  the  regularly 
attainable  feedback  gains  that  are  used  during  flight  investigations  involving 
the  AF/CAL  T-33  and  AF /CAL  C-131  (TIFS)  airplanes.  The  gains  listed  in 
the  table  are  not  necessarily  the  maximum  that  can  be  achieved;  these  gains 
are  regularly  and  easily  obtained  without  special  provision  for  sensor  and 
amplifier  noise,  servo  bandwidth,  structural  dynamics  and  other  corrupting 
or  limiting  influences.  The  gain  values  represent  day-to-day  state  of  the 
science  and  are  conservatively  estimated.  If  differences  exist  in  the  maxi* 
mally  allowed  feedback  among  aircraft,  the  more  conservative  number  is 
always  chosen.  It  is  nevertheless  important  to  note  that  large  differences 
among  individual  aircraft  do  occur  and  these  differences  are  due  mainly  to 
variations  in  structural  flexibility.  Because  a  CCV  may  be  highly  flexible, 
die  gains  listed  below  may  even  be  too  optimistic  without  specific  structural 
mode  control  provisions.  If  feedback  is  provided  to  augment  the  fundamental, 
rigid  body  dynamical  behavior  and  the  attainment  of  this  feedback  requires 
extensive  structural  mode  control,  then  the  problems  of  reliability  and  air¬ 
craft  parameter  identification  ar'  doubly  critical.  Advanced,  accurate 
methods  of  vehicle  parameter  identification  are  only  now  being  developed 
(Reference  1 1)  and  it  will  likely  be  another  decade  or  more  before  adequate 
tools  are  available  to  identify  the  structural  mode  and  flutter  parameters  of 
an  airplane,  in  addition  to  the  rigid  body  stability  derivatives. 

Two  general  comments  can  be  made,  then,  about  maximally  usable 
feedback  gains.  First,  as  mentioned  abc/e,  the  smaller  and  more  rigid  the 
airframe,  the  higher  the  attainable  feedback  gains.  Second,  the  lowe  '  the 
degree  of  the  derivatives  of  the  state  variable  used  for  feedback,  the  higher 
the  feedback  gain  that  can  be  attained.  This  limitation  is  attributable  as 
much  to  the  sensor  characteristics  as  to  the  structural  flexibility.  Pitch 
accelerometers  generate  more  noise  than  pitch  rate  gyros  which  in  turn  tend 
to  be  less  noise -free  than  attitude  gyros.  Thv,re  are  exceptions,  particularly 
when  the  good  and  bad  features  of  angle  of  attack  vanes  and  accelerometers 
are  evaluated  in  their  ability  to  alter  the  short  period  natural  frequency  of 
the  airplane  and  compared  to  the  noise  output  level  of  these  sensors,  but  in 
general,  the  state  derivative  rule  holds. 
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As  discussed  in  a  previous  section,  the  control  system  of  the  T-33 
CCV  airplane  requires  mainly  pitch  rate  and  angle  of  attack  (or  normal  accel¬ 
eration)  feedback  in  the  longitudinal  plane  of  motion  with  very  little  speed  or 
attitude  augmentation.  There  is  no  reason  to  believe  that  other  vehicles 
would  be  significantly  different  unles,s  they  possess  highly  objectionable  speed 
stability  or  phugoid  characteristics,  so  the  major  effort  was  placed  upon 
estimation  of  the  pitch  rate  and  angle  of  attack  gains,  with  little  emphasis  on 
attitude  and  none  at  all  on  velocity.  Table  X  reflects  this  emphasis. 


5.  2  CCV  BARE  AIRFRAME  STABILITY  DERIVATIVES 

Using  the  attainable  feedback  gains  given  in  Table  X,  it  is  a  rela¬ 
tively  straightforward  computation  to  estimate  the  minimum  dimensional 
stability  derivatives  that  must  inherently  be  possessed  by  the  bare  airframe 
of  a  Control  Configured  Vehicle,  assuming  that  sufficient  surface  effective¬ 
ness  and  power  is  available  to  augment  the  derivatives  if  they  are  found  to 
be  acceptable. 


The  equations  of  motion  of  the  bare  airframe  CCV  and  the  flying 
qualities  model  can  be  expressed  as 

i  =.  (64) 

(^>5) 

The  feedback  control  law  is  of  the  form 

(66) 

so  a  fully  augmented  aircraft  can  be  described  by  the  equation 

%=i{f-GK)Y.+G„u  (67) 

If  the  augmented  CCV  is  to  fly  as  the  flying  qualities  model  flies, 

then 

(68) 

ar-,'  the  matrix  of  dimensional  stability  derivatives  of  the  bare  airframe  are 
re  ,•  cted  to  the  range 

^bare  CCV  ""  (69) 
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If  the  rn'itrix  of  stability  derivatives  of  the  model,  Table  III  and  the 
matrix  of  gains  (assuming  three  controllers)  were  substituted  into  Equation 
(69)>  the  longitudinal  requirements  become: 


toartCCN 


‘■.01424- 

0 

.OOH 

-.000157 


-  32.17  0 

0  i 

0  -1.085 

0  1 


5.70 

<i,1 

%a 

0 

0 

0 

0 

-28.4 

f 

U 

Sat 

Saa 

-2.115 

9ii 

Sit 

Sta 

iS 

±4 

42 


i  J 
i3 
4Z 


±  S 
45 
44 


-,01424-4  — 
0 

.0037  4  — 
-‘.0001574- 


-32.174  59,,  4  4^2,  4.  2j,5 
0 

*  ^9,1  ^  ^9iz  -fc 
i  ^941  ^  ^94Z  ^  ^^43 


^%^^9n^29,s 

1.0 

^  1^^931-^  ^95s 

t3g^,J:3q^^4Z94.y 


S.  73  5g„  jt  Sg,g  ±  ^9t3 

0 

'ZB.4*5g^,±Sg^^^  5g„y 

-2.215 

(70)  ~ 


where  the  symbol  — ■  means  that  the  gain  is  not  considered  and  has  been 
shown  to  be  negligible. 


It  becomes  clear  then,  that  surface  effectiveness  plays  a  very  strong 
role  in  the  augmentation  possible.  If  this  aircraft  possessed  the  surface 
effectiveness  of  the  elevator,  flap  and  ailerons  of  the  T-33,  the  CCV  bare 
airframe  requirements  become 
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which  could  quickly  lead  one  to  the  conclusion  that  the  bare  airframe  stability 
derivatives  are  meaningless  to  a  CCV  if  augmentation  is  to  be  fully  realized 
and  if  the  control  surfaces  have  sufficient  effectiveness  to  do  the  job. 

If  sufficient  control  effectiveness  is  not  provided,  then  it  will  be  dif¬ 
ficult  to  realize  CCV  objectives.  Consider,  for  instance,  the  T-33  airplane 
with  shortened  tail  and  reduced  horizontal  tail  surface  =  *1'^  =  0.  5.  The 
matrices  of  stability  and  control  derivatives  become: 


If  the  matrix  of  Equation  (73*  is  compared  to  that  of  the  model,  it 
is  seen  that  two  very  important  terms  in  the  matrix  -^3^  =  and 

■P34  =  ,  which  are  used  to  approximate  short  period  damping 

and  cD^  ,  are  just  barely  attainable.  If  only  elevator  is  used  for  augmen¬ 
tation,  then  the  stability  derivatives  cannot  be  independently  altered.  Once  ^ 

a  value  of  feedback  gain  is  selected,  each  column  of  the  F  matrix  is  fixed,  I 

or  another  way  to  express  this  constraint  is  that  since  the  Ist,  3rd,  and  j 

4th  rows  of  the  F  matrix  contain  the  coefficients  of  a  separate  degree  of  ; 

freedom  of  motion  of  the  vehicle,  ? 

The  flap  and  the  aileron  of  the  T  -33  in  this  example  can  just  barely  ; 

be  considered  independent  control  devices,  because  they  produce  pitching  1 

moments,  Z  and  X  forces  almost,  but  not  quite,  proportional  to  each  I 

other.  Therefore,  normally  the  two  surfaces  would  not  be  considered  | 

efficient  for  dynamic  augmentation  purposes,  but  they  do  provide  a  good  i 

measure  of  controllability  of  the  wing  lift  distribution  and  therefore  wing  j 

root  bending  moment  control. 


5.  3  CONTROL  (FEEDFORWARD)  AUGMENTATION 

Feedforward,  or  control  augmentation  is  less  restrictive  than  feed¬ 
back.  If  the  Control  Configured  Vehicle  does  not  exhibit  the  required  con¬ 
trol  effectiveness  to  produce  good  handling  qualities,  then  other  surfaces  can 
be  used  to  augment  the  vehicle. 

The  matrix  of  control  effectiveness  terms  of  the  flying  qualities 
model  is  given  by 
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using  just  the  elevator  to  produce  the  required  control  forces  and  moments. 
The  T-33  with  reduced  tail  length  and  area  but  using  the  flap  and  ailerons  is 
given  by 


The  effectiveness  of  the  other  surfaces,  i.e.,  the  flap  and  ailerons, 
can  be  used  to  augment  the  elevator  as  shown  in  the  following  sketch: 
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Figure  13  FEEDFORWARD  OR  COMMAND  AUGMENTATION 
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The  gains  are  computed  from 


or 

Id,  =  3.7/  ,  /d^  =  -0.895' ,  =0.995- 

and  the  system  can  be  mechanized  without  difficulty. 

The  investigations  and  results  of  this  section  have  shown  that  sta¬ 
bility  and  flying  qualities  should  have  very  little  influence  on  the  geometry  of 
a  Control  Configured  Vehicle.  If  sufficient  control  power  is  available,  the 
vehicle  can  have  almost  any  shape  and  stability  augmentation,  within  the  pre¬ 
sent  state  of  the  art,  can  alter  the  flying  qualities  to  the  desirable  Level  1 
behavior. 


5.4  WEIGHT  CONSIDERATIONS 

The  T-33  configuration  with  d/t  =  Idi  =  .5  is  about  the  maximum  tail 
reduction  possible  for  the  CCV  to  still  have  the  ability  to  obtain  a  Level  1 
airplane  with  a  realizable  feedback  system.  This  corresponds  to  a  change 
in  static  margin  of  about  -15%.  The  entire  aft  section  of  the  fuselage  of  the 
T-33,  which  includes  the  tail  assembly  and  exhaust  extension  past  the  engine, 
is  estimated  by  CAL  personnel  to  weigh  approximately  700  pounds.  With 
the  .5  dfi  and  this  can  be  reduced  to  350  pounds.  Further  weight  savings 
from  a  lighter  wing  structure  for  reduced  wing  loads  due  to  the  MLC  system 
of  wing  flaps  is  conservatively  estimated  at  150  pounds.  This  500  pounds 
is  5%  of  the  normal  T-33  dry  weight. 

This  estimate,  however  does  not  take  into  account  the  weight  of  the 
reliable  control  system  that  would  have  to  be  added  to  the  existing  airplane 
to  allow  the  vehicle  to  fly  with  the  geometrical  configuration  assumed  in  this 
report.  Overall  it  is  believed  that  the  total  weight  reduction  along  with  the 
reduced  wing  root  bending  moment  and  drag  will  significantly  improve  the 
performance  or  payload  capacity  of  the  T-33. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


6.  1  CONCLUSIONS 

The  purpose  of  the  research  reported  in  this  document  is  to  investi¬ 
gate  CCV  design  and  control  system  concepts  in  a  general  way  and  to  apply 
these  concepts  to  a  T-33  airplane  in  as  realistic  a  way  as  possible  within  the 
limited  scope  of  the  program.  The  study  accomplished  this  goal  in  a  positive 
sense  and  the  following  conclusions  are  drawn  from  the  results.  These  con¬ 
clusions  are  tentative;  not  all  of  the  aspects  of  the  problem  we  e  considered 
and  further  effort  may  modify  some  of  the  results  but  probably  not  significantly 
alter  the  fundamental  principles.  The  major  conclusions  are  summarized 
below; 


2. 


3. 


4. 


5. 


Relaxed  static  stability,  maneuver  load  control,  and  good 
flying  qualities  can  be  made  to  be  compatible  if  adequate 
numbers  of  independent  force /moment  producing  devices  with 
adequate  effectiveness  and  power  are  provided. 

Because  the  geometry  and  surface  configurations  are  generally 
fixed,  the  application  of  CCV  concepts  after  the  fact,  i.e., 
on  a  presently  existing  airplane,  will  yield  only  limited 
success.  Control  Configured  Vehicle  concepts,  to  be  most 
effective,  should  be  incorporated  into  the  preliminary  design 
stages  of  a  new  airplane. 

The  T-33  configuration  with  ^  =  .  5  is  about  the 

maximum  reduction  possible.  This  corresponds  to  a  change 
in  static  margin  of  about  -15%.  Total  structural  weight 
reductions  may  amount  to  about  500  pounds  or  approximately 
5%  of  the  normal  T-33  dry  weight.  This  can  be  interpreted 
to  mean  a  10%  increase  in  fuel  capacity. 

Flying  qualities  have  a  significant  effect  on  the  control 
system  configuration.  Flying  qualities  are  very  flexibly 
or  broadly  defined  and  can  be  selected  to  benefit  the  maneu¬ 
ver  load  control  objectives  of  the  CCV,  Flying  qualities 
requirements  can  be  chosen  to  restrict  or  enhance  the 
application  of  maneuver  load  control. 

The  present  state  of  the  art  of  feedback  control  allows  for 
augmentation  of  an  extremely  wide  variety  of  bare  airframe 
characteristics  and  therefore,  geometrical  shapes  of  the 
airframes.  Stability  constraints,  such  as  static  margins, 
have  little  or  no  importance  if  sufficient  control  effectiveness 
and  power  are  available  to  provide  for  good  flying  qualities 
and  maneuverability. 
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6.  A  more  comprehensive  effort,  considering  many  aspects 
not  included  in  this  study  would  be  necessary  to  optimize 
the  results,  but  feasibility  has  been  demonstrated. 


6. 2  RECOMMENDATIONS 

This  study  has  only  investigated  a  few  CCV  concepts;  1)  reduced 
static  stability  through  reduced  tail  area  and  tail  length,  2)  maneuver  load 
control  in  a  pullup  through  the  addition  of  inboard  and  outboard  direct  lift 
flaps  on  the  wing,  and  3)  the  constraint  of  flying  qualities  requirements. 

The  feasibility  of  this  type  of  vehifcle  has  been  demonstrated.  However,  in 
the  early  design  of  any  new  CCV,  a  more  extensive  study  should  be  carried 
out.  With  a  precise  knowledge  of  mission  and  performance  criteria,  more 
specific  candidate  controllers  should  be  investigated.  Canards,  ventral  fins, 
split  wing  flaps,  and  various  wing  positions  should  be  evaluated  to  optimize 
the  design  to  achieve  desired  characteristics.  Torsional  wing  bending 
moments  must  be  evaluated,  as  this  may  become  a  problem  with  the  addi¬ 
tional  wing  flaps. 

The  concepts  of  relaxed  static  stability  and  maneuver  load  control 
represent  evolutionary,  rather  than  revolutionary  advances  in  aircraft  sta¬ 
bility  and  flight  control  practice.  Feedback  to  augment  damping  is  already 
in  full  operational  use,  so  additional  feedback  and  command  augmentation  to 
improve  static  stability  is  only  a  step  beyond  present  design  procedures. 

Wing  surfaces  designed  to  alter  the  magnitude  and  symmetry  of  the  lift  along 
the  wing  have  been  in  use  for  fifty  years  or  more.  The  use  of  flaps,  ailerons 
or  other  wing  surfaces  to  alter  the  lift  distribution  on  the  wing  during  either 
trimmed  or  transient  flight  is  als''  logical  extension  of  present  practices, 
so  maneuver  load  control  is  also  .  '  iible. 

Feasibility  only  has  been  demonstrated  by  the  results  presented  in  this 
report.  A  more  comprehensive  study  and  simulation  program  is  needed 
before  actual  mechanization  and  flight  testing  could  be  undertaken.  In  general, 
control  surfaces  are  sized  and  located  on  aircraft  only  after  extensive  anal¬ 
ysis  and  model  testing  has  been  done.  To  be  most  effective,  CCV  concepts 
should  be  included  in  the  preliminary  design  stage  of  an  airplane.  The 
application  of  CCV  concepts  after  the  fact  of  the  airplane  design  will  likely 
be  not  as  effective  and  the  modifications  will  probably  be  costly. 

Flying  qualities  requirements  will  play  an  important  role  in  the 
establishment  of  CCV  airframe  designs  and  augmentation  configurations.  Two 
flying  qvialities  models  were  used  in  the  present  study.  The  flying  qualities 
parameter  rt/cL  was  shown  to  have  a  strong  effect  on  the  use  of  flap  and  ele¬ 
vator  surfaces  to  simultaneously  obtain  lower  wing  root  headings  and  good 
flying  qualities.  In  addition,  relaxed  static  stability  of  the  bare  airframe  will 
require  relatively  large  surface  deflections  and  surface  rates  to  artificially 
produce  the  stability  characteristics  dem.;nded  by  flying  qualities.  However, 
hying  qualities  requirements  are  broad,  and  it  appears  possible  to  be  able  to 


satisfy  both  flying  qualities  requirements  and  make  optimum  use  of  CCV 
concepts  at  the  same  time.  The  mutual  overlap  of  these  requirements  should 
be  carefully  defined. 

It  would  be  very  important  to  consider  the  effects  of  geometrical 
alterations  throughout  the  entire  flight  envelope  of  operation  of  the  vehicle. 

It  is  one  thing  to  design  a  flight  control  system  that  will  give  good  flying 
qualities  at  a  single  flight  condition  but  an  entirely  different  matter  to  design 
a  simple  system  for  the  entire  range  of  operation  of  the  vehicle.  The  bare 
airframe  dynamics  and  the  optimum  flying  qualities  model  changes  as  a 
ftmction  of  flight  condition.  Although  all  the  elements  that  make  up  a  Level  1 
airplane  are  broad  in  range  at  one  flight  condition,  a  minimum  complexity 
control  system  that  satisfies  all  elements  for  a  Level  1  airplane  at  all  flight 
conditions  represents  a  formidable  challenge  to  the  Control  Configured 
Vehicle  designer. 
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APPENDIX  I 


LINEAR  OPTIMAL  CONTROL 


The  linear  optimal  control  problem  treated  in  Section  III  of  this 
report  is  a  variation  of  the  general  problem  using  the  performance  index 


(I-l) 


subject  to  the  constraint  of  the  differential  equation^  ‘'A  motion 


-  F»i  +  Gu< 


(1-2) 


y  =  Mid 


The  solution  to  this  problem  requires  that  the  Euler- Lagrange 
equations  be  satisfied 


i(i)  _  ■  F  -G1Z~^G^'\  r  %ii) 

Mi)  X(t) 


Mo)  * 

Mo) 


subject  to  the  boundary  conditions  on  the  state  vector  X  (0  )  =  and  the 
Lagrange  multiplier  X,  (0)  =  X.  .  The  basic  problem  is  to  determine  the 
boundary  condition  X  (0)  as  a  tunction  of  the  state  vector,  thereby  eliminating 
the  two-point  boundary  value  aspects  of  the  problem  which  will  then  yield 
a  closed  form  solution. 


It  has  been  shown  by  R.  E.  Kalman  (Reference  7)  and  others  that 
X  (0)  and  X  (0)  are  related  by  the  equation 


a(o)  -  y(oo)  x(o) 


(1-5) 


where  T  (oo)  is  the  positive  definite  symmetrical  solution  to  the  matrix 
Riccati  equation 


TF  f  r^v  -  pQie'V'p  h 


(1-6) 


It  has  been  shown  (Reference  10,  for  instance)  that  the  eigenvalues 
of  the  Hamiltonian  system  of  Equation  1-4  consist  of  the  eigenvalues  of  the 
stable  optimal  closed-loop  system  with  negative  real  parts  A  and  the 
eigenvalues  of  the  unstable  "adjoint"  system  -  A  with  positive  real  parts. 


If  we  transfer  Equation  1-4  into  the  diagonal  canonical  form  using  a 
linear  transformation 


Preceding  page  blank 
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. .  ,  ,.  ,  ,  .  . . ....  .  . 


ff  > 


from  which  we  can  obtain 
At..  _  At 


K(i)  =  t;,  e '^1^''^  ■"  T2 


(1-13) 


The  optimal  solution  for  ^(f)  *  which  is  stable^  cannot  contain 
terms  in  e'^'*'-  Therefore,  £^o)  and  A. (o)  must  be  related  by  the 
expression. 


A(o)  =  -  N"'  Mi:(o) 


(1-14) 


Substituting  this  expression  for  A,  (^)  1*^  Equation  1-13  yields 


lid-)  =  r„e^^  lk-lh'^m)  liio) 


(1-15) 


From  the  identity 


'T. 

’"a" 

■  K 

L' 

'  IC 

L' 

■t;. 

'  I 

“1 

.  M 

N  J 

.  M 

N. 

-Tz, 

T 

zz 

•• 

.  0 

iJ 

(1-16) 


We  can  obtain,  among  others,  the  following  relationships 


M-  -NT^,V' 

(a) 

•k:  - 

“LT  T  ■’ 

^  «  IZ 

(c) 

(b) 

L  = 

1 

1 

(d) 

(1-17) 


Substituting  these  expressions  for  fC  ,  L  ,  M  ,  and  N  in  Equation  1-15  yields 

=  1;,  ,1-18) 


Comparing  Equation  1-5  and  1-14  we  find  that  the  steady-state  solution 
to  the  Riccati  Equation  is  given  by 


T(cc)  =  -N'^M 
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(1-19) 


_ _  _ _ _ 


A  similar  development  involving  the  upper -right  partitioned  matrix 
of  Equation  1-22,  i.e.,  that 


/r/=-r  -  i'  h'qht^^  - 


(1-26) 


will  yield  another  result.  Substituting  ^  ^  ^  ^  post  multiplying 

the  entire  equation  by  -/,  *  and  pre  multiplying  by  '*’Xz  ‘^  yields 

(1-27) 


y(u>)  =  T  r  ^  T  r 

'  ;t/  //  zt  a 


(1-28) 


Therefore,  once  we  know  the  eigenvectors  of  the  Hamiltonian  system 
of  Equation  1-4  we  can  directly  compute  the  optimal  feedback  control  law. 

The  technique  shown  in  this  section  is  independent  of  the  order  of  the 
system  and  the  number  of  controllers  of  the  dynamic  system.  Equations 
1-18  and  1-20  show  that  the  feedback  gains  and  the  regulator  transient  response 
are  directly  related  to  the  eigenvectors  of  the  system.  Those  eigenvectors 
are  a  function  of  the  weighting  matrices  Q  and  R  in  the  performance  index. 

Example 

Consider  the  two  controller,  second  order  system  described  by  the 
equations 


^  ""1  f  "^'1  ^  f'  '’If 

^  I  /  '3}  I  %  \  *  lo  i  \  I  U,- 


(1-29) 


It  is  desired  to  find  the  optimal  control  law  that  will  satisfy  the  performance 
index 


(1-30) 


where  ^  X  "  X  ^  ^ 


The  Hamiltonian  system  for  this  example  is 

r-z  z  -i  0  r  X.,  - 

‘^z  f  -3  0  -1 

j,  "  1  0  *2  1  X, 

-A,  J  L  0  i  2  ^3  J  L  A, « 


(1-31) 


whose  characteristic  polynomial  is  given  by 


A(s)a{s)  -  *  S.56  5  +  S.S9)  ($*■  -^.SS  %  +-  5',S’9) 

=  (S  -  4.116')  (5  i 


(1-32) 


The  left  half  plane  roots,  i.e.,  the  eigenvalues  of  the  closed-loop 
optimal  system,  are  given  by 

S,  ,  ”  -  ,  -4.  Il< 

and  it  is  necessary  to  find  the  two  eigenvectors  of  Equation  1-31  associated 
with  these  two  eigenvalues 

The  eigenvector  transformation  is  found  to  be 


r  = 


.  Ai>^  -•  ^3^  I  ♦  * 

.it?"  ".747  '  ~ 

'27/  ',0^3  I  '  ‘ 


"7/  I  "Tt 
"T-f  I  "7^ 


(1-33) 


where  the  two  blank  columns  represent  the  eigenvectors  associated  with  the 
eigenvalues  ^  =  ■i’f.44,  -f  4.  From  Equation  1-25  ,  we  have 


which,  for  our  example  becomes 


( 

t 

\ 

■  e  0  ^ 

.4.IIS  t  1 

.  0  t 

- 

•IL,  (t)  -  .  3Z2 

^  (t)  ^  (.3Zb  %,(0) 
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•7  2 


continued 


t _ , 

INTEGRATE  FOR  ONE  TIME 
POINT  AND  CALCULATE 
NEW  STATE  VECTOR: 

X,  AT  PRESENT  TIME  POINT 


YES 


PRINT  OUT  FINAL  VALUES 
AND  TIME  HISTORIES 


INPUT 


All  input  is  read  in  NAMELIST  form: 

u,(^3)  =  initial  input  steps  of  Sg  ,  5*  ,  Sq 

-  model  F  matrix 
=  weighting  matrix  on  states 
XL(3)  =  initial  values  on  2»  3 

TITLE  =  (any  80  characters) 

=  constants  for  F,  G,  1 

matrices: 


F 


Z(l)l<:^fZ(Z)  -32.11  0  Z(3)Zf,  +  Z(4) 

0  0  1  0 


Z(11)  0 


z(iz)i(:j^>z^tz(i3) 


z(r^) 

0 

0 

9^(17) 

i(zo) 

^(21) 
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0)1 

NTP 

V 

T 

CA 

■REDUCE. 

LOOP 

NCVCLE 

yjSTOP 


IPRINT 

KA 

KL 


=  ^  time  increment 
=  number  of  time  points 
=  weight  on  drag 
=  weight  on  WRBM 

=  initial  constant  for  cC>i  (defaulted  to  1.) 

=  reduction  multiple  for  update  (defaulted  to  .  1) 

=  max  number  of  loops  in  updates  (defaulted  to  10) 
=  max  number  of  conjugate  gradient  iterations 
(defaulted  to  10) 

=  J  convergence  test:  stops  if 
(defaulted  to  .  005) 

=  to  print  intermediate  interations  set  to  1 
if  not  set  to  0.  (defaulted  to  1) 

=  initial  tail  area  ratio 
=  initial  tail  lengtlx  ratio 


OUTPUT 


The  output  form  is: 

(  Title  Date 

Final  values  wanted  for  states:  ^,04, 
Model  1  matrix 
Q 


input 

data 


Zi  Drag,  WRBM 


•  ^2  *  ^3 
F  matrix 

Linear  search  values  for 


OC/ 
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I 

I 

t 

•I 


fy.^j  ■  '':^'5,'--  '  ?i?s 


CONJUC.ATF  GRAniFNT  MPTHOO  FQK  f.CV  tTH  FIND  TAIL  ARE  A  t  LFMGTH,  AND  CDMSTAK'T 
control  inputs  DE.ni,DO  FOR  A  4G  PULL  UP 

DIMENSION  X(«,l) ,XOUT(4,Sl »,F(4.«»  tGI A ,3 ) ,U I  3 » ,L I  4, 4 ) t C «  3) t Q( 4 ) , 

♦XL(  3)  fTIILFIPO)  tFFOM(4,4)  ,Fr.nM|4t3>,Gr.nM(  3,3)  ,f)FFA(4,4),n«=FL(4,4) , 
*DFGA(4,3)fDrr,L(4,3)tnC,GA(  l,3),nGGL(  l,3),n2(3,n,M?(3f  1  ),FMLC4,4), 
♦SRI4f3) fPHIl(4,4)  ,PHl2<4,4»t2(40),TIME(5n ,nF A { 4, 4 ) ,nFL( 4 ,4 » , 

♦  DGAtAtl)  ,TM2M2(3t3),TEMPH4,4)«TCMP?(4,4),TFMP3(4,4),FMLT0(4,4) 
*,nFAT0(4,4) ,nFLTg(4,4»,SKA(4,l»,SKL(4,l) ,PHI2G(4,3»,MHI?GU(4tH  f 
♦DRAGISl)  tWP9M(bn,XTIlt4).,TFMP4(4,4),SKAT(  I,4),SKLT(  1,4)  ,SRTn,4)  , 

♦  A(B)  ,APPFV(ft)  ,ULAST(3) ,XLLAST(3) , DEL JIB) 

REAL  XA,KL,L,M1 ,M2,KALAST,KLLAST 
RFAL*8  TOATE 

NAMFL!ST/INPUT/U,L,C,0,XL,TITLF,Z,nT,NTP,V»TtCA,RFDUCF,Lnr:P,NCYr.i  F 
♦tXJSTOP, IPRINT,KA,KL 
CALL  CLFAPICA.OELJIfi) ) 

CA=l. 

XJSTnP=.005 

REnuCF=.l 

Lpnp=in  - - - — ”1®^ 

NCrCL-lO 

1PRINT=1 

CALL  OATFITOATF) 

I  REAOIb, INPUT, END=99D9) 

WRITE(6,20)TITLE ,rPATF,(C( I),I=1,3),I(L(I, J),J=l,4), t-l,4) 

20  F0RMAT(1H1,2X,20A4,10X,A8,//,?X, ‘“OnFl  0 , ALPHA , Al  OHA-DOT  AT  PJNAL 

♦  TIME  ARF:  •  ,3E17.A,//,2X, ‘MnOFL  L  MATRIX  IS : • , /, 1 4E16 .4  )  ) 

HRITF(6,2I)  (0( 1), 1*1,4), V,T 

21  FORMAT!//, 2X,‘0  =  •  ,4F17.6,/,?X*V  *  •  ,Fl  7.6, /,2X,  •  T  =  SFW.A) 

IST0P*0 

ICYCLF=0 

SET  UP  VARIABLES  WHICH  ARF  NOT  FUNCTIONS  OF  KA,KL 
FI  1, 2)*- 32. 17 
FI2,3)=1. 

FI4,l)=2Ill ) 

F(4,^'*l . 

GIl,^  .=ZIIS) 

Gtl,3)=Zn6) 

GI3,2)=Zn8) 

GI3,3)=ZI!9) 

GI4,2)*Z(21) 

GI4,3)=ZI?2) 

02I2,l)*ZI26) 

n2{3,l)=ZI27) 

M1»ZI2R) 

D01A=ZI24) 

nD2A=ZI?S) 

M2I2,l)  =  7.I29) 

M2I3,l)*ZI30) 

0FAIl,l»^ZI2) 

DFAI l,4)*2I4) 

DGAI 1,1)*7I 14) 

DGAI4,1 )*7120) 

DO  30  1=1,3 

30  TEMPII 1 ,I )*M2I I,) )*T 

CALL  MATMPYIM2,TFMP1,TM2M7,3,1,3,3,4,3) 

35  ICyCI.F*ICVCLF+l 

36  CONTINUF 
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C  SET  UP  TERMS  WHICH  nFPFNn  UN  KA,KL  NFEPFD  FOR  J  CALCIJLATfONS 
F(l,l)=Z(l)^?(2>*KA 
F(l,4)=Z(3»*Z(A>*KA 
F(3tl)=Z(5)+Z(f.)*KA*K«. 

FO,3)=ZI7)^Z(8)*KA*KL 

F(3t4»=Z(«))tZ(lO)*KA*KL 

F<4,4)=Z(l2)*Z<13i*KA*XL 

R(l,l>=ZIl4)*KA 

G(3tn=Zn7)*KA*KL 

G(4,1)=Z120»*KA 

ni=Z(23)+Z(24)*KA 

n2(  l,n=Z(?5)«‘KA 

CALL  MATFXPt4,nT.F,PHtl,PHI2#TFMPl,TFHP2,4»4,4,4,4,LT» 

IFILT.FO.OIGO  TO  40 

WPlTF(4,37)Lr,0r 

37  FnRMAT()H0,3X,'MATFXP  FATIEO,  LT  =  M2t*  OT  =  SFfl.b) 

3fl  WR1TE<6»39)( (F( If J)t J=lt4),I=l,4) 

3R  FORMATIlHOf 3X,*F  MATRI X  * t / 1 <4F16.6) ) 

GO  TO  I 

40  no  41  I=?vNTP 

41  T!M?( I )=TIMF( I-l J+OT 

CALL  MATMPY(PHT2fG»PHI2Gf4,4t3,4,4,4) 

CALL  HATMPYI  PHIPGfU,  PHIpGUt 4,3*1  t4t.3,4) 

CALL  3ATA00(F,L,PML,4,4,4,4,4,n 
C  SFT  UP  haTPICES  WHICH  OEPFNO  ON  KA,KL  ,  FOR  J  CALCULATIONS 
00  45  1=1,3 

45  TEMPl(l,t)=0?n,l  J*V 

tall  MATMPY(n2,TFMPl,TFMP2,3,l,3,3,4,4) 

CALL  HATA00(TeMP2,TM2M?,TEMPl,3, 3, 4, 3, 4,0) 
no  46  1=1 ,4 
no  46  J=l,3 

46  TEMP2(  J,n=G(  l,J)*on  1 

CALL  MATMPY(TFMP2,G,TEMP3,3,4,3,4,4,4) 

CALL  MATAOni TEMP3, TFMPl ,CGOM, 3,3 ,4,4, 3, 0 ) 

00  47  1=1,4 
00  47  J=l,4 

47  FHLTQI f ,J)=FML( 

CALL  MATMPY(FMLTg,FML,PFnH,4,4,4,4,4,4) 
FFi)M(4,4)=FFPM(4,4)  >V*ni**2+T*Ml»4? 

CALL  MATHPV(FHLT0,G,PG0V,4,4,3,4,4,4) 

FG0M(4,1  )=FGl)H(4,  n+V*Dl*n2n  ,  l» 

FG0H(4,7I  =  FGnM(4,2M\f*ni*n?(2,U»I*Ml*M2(?,l  ) 

FGOMI  4,3)  =  Fr.nM(4,3)  +  V»01*02f  3, 1)  ♦T*M1*M2(  3,1) 

IF( lOELJ.FO.O)  GO  TO  60 
C  SET  UP  MATRICES  FOP  DFLJ  CALCULATIONS 
OGAI 3,’ )=Z(17I*KL 
OGL=ZI 17:*KA 
OFAI 3,1 )=2(6)*KL 
nFAI3,3)=2(6)*KL 
OFAI 3,4) =Z( 10)*KL 
nFA(4,4)=Z( 13)*5L 
OFLn,l  )=716)*KA 
0FLI3,3)=Z(R)*KA 
0FL(3,4)=2(  ini*t(A 
0FL(4,4)=Z( 13)^KA 

on  50  1=1,4 

50  TFMPK  1  ,I  )=0GAn  ,  I  )*0(I  I 

CAlL  MAYMPY(TFMPl,G,0GGA,l,4,3,4,4,l  ) 


t.''''-*'^--  '"/'y'"^'^^'  •».’s 


DO  SI  I  =  U3 

51  OGGA(l,n*DGCA(l,n*V*00?A*0?ll,l) 
no  52  1=1,3 

52  0GGL(l,ll=nGL*0(3)*G(3,I» 
on  53  1=1,4 

no  53  J=l,4 

53  OFATO(I,J»=OFAtJ,n-tO<J» 

CALL  MATMPY(nFATQ,FML,DFFA,4,4,4,4i4,4) 
OFFA(4,4)=OFFA(4,4»»V»nniA*01 
no  54  1=3,4 
no  54  J=l,4 

54  nFLTOU,U=OFL(  !,J)*0(!) 

CALL  MATMPYIOFLTQ,FML,OFFL, 4, 4, 4,4,4,41 
CALL  MATMPY(0FATQ,G,0FGA,4,4,3,4,4.4l 
CALL  MATMPY(FMLTQ,0GA,TFMPl,4,4, l,4,4,4) 
no  55  1=1,4 

55  DFGA ( I , n  =OFGA( 1,1)  +  TEMPI ( I  ,  I » 

OFGA(4,ll  =  DFGA(4,n>V*001A*02(  I,  1) ♦V*D1*D02A 

nFGA(4,2>=0FGAt4,2l+V*0niAv02l2,l> 

0FGA(4,3l=0FGA<4,3UV*nniAvr2';3,l) 

CALL  MATMPY{OFLTO,G,nFGL,4,«,Hf4,4,4) 
no  56  1=1,4 

56  nFGLH  ,  n=OFGL( ! ,  1  l*FMLTO<  1 ,3I*DGL 

60  CONTINUE 

C  INITIALIZE  FOR  TIME  =  0. 

XJ=0. 

on  61  1=1,8 

61  06LJ(n=0. 
no  62  1=1,4 
X(  l,U=0, 

SKA( 1,11=0, 

sKLn,n  =0, 
on  67  J=l,3 

62  SB( I,JI=n. 

SUM0=0. 

suMw=n. 

no  100  IT=2,NTP 
1F( inFLJ.EO.Ol  GO  TO  70 
C  CALCULATE  SENSITIVITIES 

CALL  MATMPY<nEA,X,TFMPl,4,4,l,4,4,41 
no  6S  1=1,4 

65  TEMPl(I,l)=TEMPl(l,l ) ♦OGA ( I , 1 1 *01 1 ) 

call  MATMPY I  phi  2, tempi, temp?, 4, 4, 4, 4, 4, 4) 

CALL  MATMPYIPHll,SKA,TEMPl,4,4,l,4,4,4> 

CALL  MAT AnO( TEMPI, TEMP?, SKA, 4, I, 4,4,4,01 
CALL  MAT'«PY(nEL,»,TEMPl,4,4,l,<,4,4) 

TEMPU3,n»TFMPl(  3,1  )^OGL*uai 

CALL  MATMPYI PHI 2, tempi, TEMP2, 4, 4,1,4,4,41 

CALL  MAT MPV (PHI  I, SKL, TEMP  1,4, 4, 1, 4, 4, 4) 

CALL  MATAOO(TFMP1,TEMP2,SKL,4,1,4,4,4,0) 

CALL  MATMPV(PHI1,S3, TEMPI, 4, 4, 3, 4, 4, 4 5 
CALL  MATAOO'TEMPl ,PHI?G,SR,4,3,4,4,4,Pl 
70  CALL  NATMPy(PHIl,X,TFMPl,4,4,l,4,4,4) 

CALL  MATA00(TFMP1,PHI2GU,X,4,1,4,4,4,0) 

IF( inCLJ.EO.O)  GO  TO  75 
nn  72  1*1,4 
72  XOUTn,IT)*X(I,n 

DRAG(TT»=0lPX(4,n^n?Il,l  l*UCl)»n?l2,ll*U(  2U02(  3, 1 ) *0(3 1 
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W?^V>W39f«V'-.,iSj''^p*»^^  V  ,  ,, 


SUMn»SUM(HORAr.(  f  T  > 

WRBM(  IT)  ^Ml*X(4f  I  )fH2l?f  1 1 *U( 2)*“2I  3#  1 )  ♦'!(  3) 

';UMW=SUMW^WRRM(  IT) 

75  00  75  1=1,4 

76  *Tn,I)  =  XI  1,1  ) 

C  INCREMENT  J 

CALL  MATMPy|FFnM,X,TEMPl,4,4, I, 4, 4, 4) 

CALL  MATMPV|FRDM,U,TPMP2,4,3, 1,4, 3, 4) 
no  77  1=1,4 

77  TFMP3(  1,1  }=2.*TFMP?n  II 

CALL  HAT AOni TcMP3, TFMPl ,TFMP4, 4, 1,4,4,4,01 
CALL  MATMPVl XT, rEMP4, TFMPJ, 1,4,1 ,l ,4,1 ) 

XJ=XJ*TFMPJ*Or 
IFtIOFLJ.fQ.nl  GO  TO  0'' 

c  INCREMENT  OFLJ  fpeproduced  from  Wk 

00  ao  1=1,4  I  best  available  copy^_%p 

SKAT(l,t  l=SKAC  !,n 
SKLT(l,l  )=SKL(  1,1) 
on  80  J=1 ,3 
80  SRTl J,1 )=S01 I ,J) 

CALI  MATAnO(TFMPl,TEMP2,TFMP3,4, 1,4, 4, 4,0) 

CALL  MATMPY{SKAT,T€M03,TrMnjA,l,4,l,l,4,l ) 

CALL  MATMPY(SKLr,TFHP3,TFMnjL, 1,4, 1,1,4,) ) 

CALI  MATMPYI SRT,TFMP3,TFMP4,3,4, 1,3, 4, 4) 

CALI.  MAT»«PY(nFF4,X,TFMPl  ,4,4,1,4,4,4) 

CALL  M.ATMPY(nFGA,ll,TF‘*P2,5,3,  l,A,3,4) 

CALL  “A’’  AOOl  TE^Pl  ,  TF  HP?  ,  T  FMP3 ,4,  1,4,4,4,01 
CALL  MATmpyI XT, TeHP3,TFMP, 1,4,1, 1,4,1 ) 

0FLJ(4)=nFL  Jt4)*(TFMnjA^TFMP)<-2.*OT 

CALL  MATMPYIPFFL,X,TrMPl,4,4, l,4,4,4) 

CALL  MATMPY(nFr.L,U,TFMP2,4,3, 1,4,3,4) 

CALL  •'ATA00ITFMPl,TF''P2,T»-MP3,4,  1,4,4,4,0) 

CALL  MATMPYI X T , TFMpi , TFMP , 1 ,4,1,1,4,11 

OFLJl  5)=nEL  J(  5  l^-l  TFMnjL»TFMP)*2.*l)T 
call  MATMPVC  XT  ,Fr.n|«,  TFWPl  ,  I, 4, 3,  1,4,4) 
on  85  1=1,3 

85  nFLJ15»T)=0EI.J(5+I)«'(  TF’^P^.l  I  ,  1)  ♦  TFMP  1 1  I  ,  I  )  )  *2  .  *0T 
90  CGNTINU*^ 

100  CONTINUE 
C  FINAL  SUM  UP  OF  J 

CALL  MATMPYI  CGr)M,U,  TFMPl  ,X,  3,  1,3, 3, 4) 

XJ  =  XJ*TIVFlNTP)*(Ull)*flM?l(l,l)»-U<?)*TEMl>l<2,l)*UMI<=TFMPiri,l)‘ 
XA0r)=XL(l)*lX(3,l)-Cll))P*2«-XL(?)*lX(4,l)-CI2))**2 

Anr)  =  l/  11  1  )  *xn  ,  I  )  ♦■X  (  3,  I  )^F  (4,4)  ♦Xl  4,  I  l+Gl  4, 1)  ♦Un  )♦•"»  I  4, 7  )*U1  ■* )  ♦ 

*G(4,3I*U13)-C(  3)  )‘'*2*xl(3) 

X j=y j*xAnnfAD0 
IF{ IDFLJ.CO.O)  GO  TO  190 
C  FINAL  SUM  UP  OF  OELJ 

0FLJCI)=(X13,1 )-C(l> )**2 
nFLJ(2)=(X(4,l )-C(2) )**2 

nFLJ3  =  71  11  )YXn,l)»X(3,  l)»F(4,4)*X(4,l  )*G(4,  DfUd  )>GC4,2)*U(2) 
tFGI4,3>*U13l'CI3) 

DELJl 3) =nFL J3**? 
xLAnn=o, 

XA4nD=0. 

on  105  1=1,3 

XAAI)n=XAAnDF2.*TIMK(NTO)*Ull  )*OGCAn  ,  1  )RU'  I  ) 

105  XL  ADn=XL  A0nF2  .*T  I  MF  (  NTP)  *01 1 1  *UGGI  (l,n*UtI) 
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B  -  -- . 


0ELj(4»*nFLJUI*XAAnnfXLt  n*SKA!3f  ll*2.*U(3,  n-C(l))*Xl(2»* 
*SKA<4,1  )42.*{X(4,l  J-C<7n 

0FLJ(4)  =0eLJ{4)+XI.(3  )*(2(  in*SKACl,  n*SKA{  3»  I )  ♦F  (  4, 4  )  *SK  A(  4,  1 )  )  ♦ 
♦2.*{nELJ3) 

nFLJ(4)=nFLJ(4)fXLi3)*(nFA(4,4)*XU,  lKnGA(4,  I  )*IJ(  1)  »*2.*nt-LJ3 
nCLJ(5»=0FLJ(5J*-XLAnO«-yL(  I  )*SKLt3,ll*2.*(X(3,  I  )-C <  I )  )  ♦XL  ( 2  )  * 
i  ♦SKL(4,l»*2.*(X(4,n-C(?)i 

I  f)FL  J(‘>»=nELJ(S)+XL(^)*{7(  IH*SK.L  ( 1,  l)^SKLn,i  ^♦F(4,4)*SKL^4,1)  »  * 

♦2.*nELJ3 

i  OELJ(5l=nELJ(5»fXLt3)*r)Fl  {4,4  )*X  (4,  I  ^♦2.*nFLJ3 

*  CALL  MATMPYJGGDM, U, TEMPI ,3,3, 1,3,3, 4> 

i  nn  110  1=1,3 

110  nELJ{US»=0FLJ(I^S»  +  2.«'TIHE(NTP»*TEMPllI,n 
nn  115  1=1,3 

115  nFLJ(  I  +  5)=nELJ(  I  ♦5>^XL(  1  )*S0(  3,1  J*?.*{X{?,1  )-C(  1  (7)*SR(4,  I  )* 

*2.*1X(4,  1)-C(?n 
on  120  1=1,3 

OFLJf  U5>=0FLJ(  {♦5)^XL(3I»(2(  1 1)*S3(  l,n^SR(  3,  nfF(4,4>*S«(4,I )  ♦ 
iG(4,ll l*?.*nELJ3 
120  CGNTINIJF 

WRITF(6,1?5»T  !TLr-,TDATF,  lCYaE,XJ 

125  FORMAT  ( IHl,  2X,20A  •  ,10X,  A3, //,4X,  MCYr,LF  =  ' ,  13,4X,  •  J=  SP^.h,//) 

1F(  ir.YCLE.EQ.l  1  Gil  TO  135 
C  TEST  FOR  J  r.ONVFRGFNCF 

lF(AHS(XJLAST-XJJ.LF.ARS{XJSTnP*XJlASTn  GO  Tt)  140 
1F( ICYCIF.GF.NCYCLF)  go  TO  130 
1F(  IPRlNT-EO.n  Gf’  TO  nS 
on  rn  i50 

130  IST0P=1 

lTER  =  lCYr.LF-l 
WRITF(6,131  HTER 

131  FnRMAT(lHO,?X,‘J  FAlLEn  TO  SATISFY  CONV:-";  GENf  E  TEST  AFTFp  ',13,'  I 
iTFUATIONS.  FINAL  TIME  HISTORIES 

GO  TO  145 
135  ITFRsICYCLC-l 

WRITE(6,i361ITPR 

135  format ( IH0,2X,'AFTFR  ',13,'  ITERATIONS,  T r Mf  HISTORIES  APF:',//) 

GO  TO  145 

140  1ST0P=1 
ITER=1CYCLE-1 
HRITFI S,I4] I ITFR 

141  FORMATI  lHn,?X,' J  SATISFIED  r.ONVERGFNCF  TEST  AFTER  ',13,'  lTEBA:Io,< 
$S.  FINAL  TIMF  histories  ARE:', //I 

145  WRITFIf.,146.) 

146  FORMAT! 3 X,' TIME' ,7X, 'OFLTA  V ,10X,*TMFTA' , I 3X, ‘ O' , 1 3X, ' Al PM A*, 12*, 
♦'DRAG' ,12X,'WRBM' ,/J 

HRl  TEI  6,  14  7)  I  TIMFI  I  ),  (  XOOTU,!!  ,  J  =  1  ,4J  ,ORAGI  I  1  ,WRBM  (  1)  ,  I  =  l,NTP  » 

147  F0RMATny,F6.2,6ri6,6) 

mR  I  TEI  6, 1481  SUMD,SUMW,KA,KL,  lUI  11,1=1,31,1X1. 111,1  =  1,3) 

I4fl  FORMATI IHO,'  INTEGRAL  OF  DRAG  =  ',£13.6,/,'  INTFGRAl  OF  WRRM  =  1, 
*E13.6,//'  KA  =  ',F13.6,/,'  KL  =  •,F13.6,//,'  OF  =  ',613.6,/,' 

♦  01  =  ',613.6,/,'  00  =  ',613.6,//,'  LAMHOA  =  ',3F17.M 

WPITFI6,3y»  I  I  FI  I, J) ,J= 1,4), I =1,4) 

IFINCYCLc.EO.I )  GO  TO  I 
IF! ISTOP.FO.l )  GO  TO  1 
150  nNTINUE 

C  CONJUGATE  GRADIENT  UPDATF  I  OOP 
Br)=BN 
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■•OSiV.v 


littflit'i'i-tn'i 


irT'JiAA 


nr.  ,55  i^itB 
155  BN=PN*nFLJ( ! )*«2 

IF{  ICVCLE.eO.n  r,0  to  I6l 

RErA=RN/Bn 

DfJ  160  !=l,8 

160  A(  I  )=OFLJJ  I)«-RFTA*APHEV<I  ) 

GO  TO  165 

161  no  16?  f=l,8 
16?  A( I )=OELJ( I  I 

C  ONF  OIMFNSIONAL  SEARCH  LOOP  FOR  ALPHA 
165  !DELJ=0 

ALPHA=CA/S0RT1BN) 

XJLAST=XJ 

kalast*ka 

KLLAST=Kl 
on  165  T=l,8 

168  APRFVn)=A(Il 

nn  169  1^1,3 

ULASrC I l=U( I ) 

169  XLLASTl I 1=XL( 1 > 

170  r.OMTINUF 

on  210  iLonp=i,LnGP 
175  KAsKALAST-Al PHA*A(41 
KL=Kl.t.AST-ALPHA*A(51 
on  150  f=l|3 

U(  n=ULAST(  I  )-ALPHA*A(5*i  , 

180  XL(  n*XLLASTm-ALPHA*Am 
WPlTF(6fl55)  ILOOP, ALPHA 

155  PnRMAT(////,2X,* IN  ONE  OIH.  SEARCH  IHOP,  ! I OOO  =  SI3,'  ALPHA  =  S 
*E13.6» 
on  TO  36 

190  WRITE(6fl95l  XJ 
195  FnRHAT(4X,«J  =  •,fl3,6l 
IF(  ILOOP. EQ.l )  GO  TO  200 
IFiyj.GE.XJLOnpi  GO  TO  ?30 
GO  TO  205 

700  IFIXJLAST.lt. 0.1  GO  TO  201 

IF(XJ.LT.i.00001*XJLASTl  GO  TO  205 
GO  TO  203 

201  IF<XJ.LF..999999*XJLAST»  GO  TO  205 
?03  ALPHA=ALPHA*REnuCE 

IF! ALPHA.lt. l.F-20)  GO  TO  1 
GO  TO  175 
205  XJLnnP=XJ 

AL?HA=ALPHA*?. 

710  CONTINUF 

715  KA=KALAST-ALPHA*A(41 
KL=KLLAST-ALRHA*A(5I 
no  220  1=1,3 
Um  =ULAST(  I  )-ALPHA*A(5>l  1 
720  XL(  n  =  XLl  ASTI  I  )-ALPHA*A(  1  ) 

CO  TO  235 

?i0  AL '>HA=  I  AlPHA/?.  )*(XJ-4.*X JL00P»3.*XJL  AST  )/<?.♦  XJ-4,*XJL0nP*2.*XJ|  A 
♦  STJ 

GO  TO  215 

235  WR1TE(6,240)  ALPHA 

740  FORMAT!////, 2X,*CPMPLETF0  ONE  OIM.  SEARCH  LOOP,  ALPHA  =  SC’,  3.6) 
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